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ABSTRACT 



O ' Gravitational lensing by the dual cusp (^-3) catastrophes of the cold dark matter (CDM) 
caustic rings at cosmological distances may provide the tantalizing opportunity to detect 

O 

CDM indirectly, and discriminate between axions and weakly interacting massive particles 



(WIMPs). The infall of CDM onto isolated galaxies such as our own produces discrete 
number of flows and caustics in the halo CDM distribution. Caustics are places where the 
CDM particles are naturally focussed. In the CDM cosmology, caustics in the distribution of 



dark matter are plentiful once density perturbations enter the nonlinear regime. Our focus 



is upon the caustic rings which are closed tubes whose cross-section is an elliptic umbilic 
(-D-4) catastrophe with three dual cusps. We call this cross-section a "tricusp." Caustic 
rings are located near where the particles with the most angular momentum in a given inflow 
reach their distance of closest approach to the galactic center before going back out. We 
reparameterize the CDM flow near a caustic ring, and obtain the ring equations in space, as 
single valued functions of an angular variable parameterizing the tricusp. A caustic ring has 
a specific density profile, a specific geometry and, therefore, precisely calculable gravitational 
lensing signatures that vary on the caustic surface, depending on the location of the line of 
sight. We calculate the image magnification of point-like sources as a function of the angular 
parameter of the tricusp for the line of sights that are parallel to the galactic plane of the 
caustic ring and near tangent to the surface. The magnification monotonically increases 
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as the line of sight approaches to the cusps where it diverges in the limit of zero velocity 
dispersion. If the velocity dispersion is small but nonzero, the divergence is cut off, because 
the caustic surface gets smeared over some distance in space and the cusps are smoothed 
out. The lensing effects are no longer infinite at the cusps, but merely large. To estimate the 
magnification numerically, for the rings at cosmological distances and for the ring closest to 
us (the fifth caustic ring of the Milky Way), we choose sample points near the cusps. We find 
that the lensing effects are largest at the outer cusp which lies in (or, near) the galactic plane. 
Around the outer cusp, the surface is curved toward the side with two extra flows. We call 
such a surface concave. In the limit of zero velocity dispersion, at a sample point near the 
outer cusp, we find 37% magnification for the CDM caustic rings at cosmological distances. 
This is about 250 times stronger than the lensing effect we predicted for the concave folds 
(considering the A 2 catastrophes enveloping the galaxy), and 37 times stronger than the 
largest lensing effect we predicted for the CDM caustics (considering the convex folds of the 
rings) before. In the presence of finite velocity dispersion, the lower and upper bounds of the 
effective velocity dispersions of the axion and WIMP flows in galactic halos may be used to 
constrain the lensing effects at the smoothed cusps. For a cosmological axion caustic ring, we 
find that the magnification may range between 3% and 2800% at the outer cusp, and between 
2% and 46% at the non-planer cusps. For a cosmological WIMP caustic ring, on the other 
hand, we constrain the magnification between 3% and 28% at the outer cusp, and between 
2% and 5% at the non-planer cusps. Because the upper bounds for the magnification at the 
cusps are obtained considering the minimum primordial value of the velocity dispersions for 
the CDM candidates in space where no small scale structure has formed, and also because 
the observer's line of sight may not exactly be parallel to the galactic plane of the caustic 
ring in general, they should be regarded with precaution. The images of extended sources 
may also show distortions that can be unambiguously attributed to lensing by dark matter 
caustics. Finally, we derive the Catastrophe Function of the triaxial caustic rings. We obtain 
the flow equations as the equilibrium points of this Catastrophe Function. The analysis of 
the Stability (Hessian) Matrix show that the caustic rings are structurally stable. 



PACS numbers: 98.80.Cq, 4.62,+v 
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I. INTRODUCTION 

A large amount of astronomical evidence indicates the existence in the universe of more 
gravitationally interacting matter than luminous (baryonic) matter. Most of the matter 
in the universe is dark (23% of the composition, according to the first year data [1] of 
the Wilkinson Microwave Anisotropy Probe, is dark matter, whereas only 4.4% is ordinary 
baryonic matter). The stuff that is responsible for holding galaxies, and clusters of galaxies 
together, is a peculiar kind of matter that we neither see nor detect by any means. The 
detection dilemma of dark matter particles has been occupying physicists for many decades 
by now, and the case still remains open. There are, however, significant reasons to believe 
that the dark matter of the universe is constituted, in large fraction, by collisionless particles 
with very small primordial velocity dispersion. Such particles are called cold dark matter 
(CDM). The attribute "collisionless" indicates that the particles interact so weakly that, 
practically, they move purely under the influence of gravity. Hence, unlike ordinary baryonic 
matter, they are undissipative. The small primordial velocity dispersion allows CDM to 
clump in galactic scales. It is believed that galaxies are surrounded by CDM which keeps 
falling into the gravitational potential wells of galaxies from all directions, and forms halos 
around their baryonic disks. The leading CDM candidates are axions and weakly interacting 
massive particles (WIMPs), such as neutralinos. The information about the distribution of 
CDM in galactic halos is substantial for the detection experiments, for the understanding 
of galaxy formation and galactic structure. The CDM infall onto isolated galaxies, such as 
the Milky Way, produces a discrete number of flows (the minimum number of CDM flows 
on the Earth is of order 100 [2]) and associated caustics throughout the halos of galaxies. 
Caustics are places in the halos where the CDM particles are naturally focussed and, hence, 
the density is very large. 

Zel'dovich [3] emphasized the importance of caustics in large scale structure formation, 
and suggested using the name "pancakes" for them. The reason why galaxies tend to lie on 
surfaces [4] is that the 3D sheet on which the dark matter particles and baryons lie in phase 
space acquires folds on very large scales, producing caustics called Zel'dovich pancakes. 

Sikivie [5, 6] derived the minimal CDM caustic structure that must occur in galactic 
halos. There are two types of caustics in the halos of galaxies: inner and outer. The outer 
caustics are simple fold (A 2 ) catastrophes located on topological spheres surrounding the 
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galaxy. Outer caustics occur where a given outflow reaches its furthest distance from the 
galactic center before falling back in. The inner caustics are rings [5, 6]. They are located 
near where the particles with the most angular momentum in a given inflow reach their 
distance of closest approach to the galactic center before going back out. A caustic ring 
is more precisely a closed tube with a special structure. Its transverse cross-section is an 
elliptic umbilic (-D_ 4 ) catastrophe which is a closed line with three dual cusps, one of which 
points away from the galactic center; see Fig. 1. We call it a "tricusp." We use the terms 




FIG. 1: Cross-section of a caustic ring in the case of axial and reflection symmetry, and p,q <C a. 



"cusp" and "dual cusp" interchangeably. In this paper, our focus is upon the caustic rings. 

The motivations to study CDM caustics originate from their analogs in Optics. Caustics 
are well known phenomena in light. They are singularities in the propagation of light. 
Examples are the bright wavy lines generated at the bottom of a swimming pool by the sun 
light, or cusp-like reflection (cardioid) that is sometimes seen in a coffee cup. Generically, 
caustics are boundaries between the two regions where the light intensity decreases very 
sharply on one side, and, typically two light beams superimpose on the other. Two conditions 
must be satisfied [7] for caustics to occur: (1) the propagation must be collisionless, (2) the 
flow must have small velocity dispersion. Propagation of light is classically collisionless, and 
the flow of light from a point source has zero velocity dispersion. Cold dark matter flow 
also has both of the properties that light has, i.e., (1) CDM particles interact only weakly, 
and (2) the primordial velocity dispersion 8v of the cold dark matter candidates is of a very 
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small order [6, 10]. For the axions, we find 
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where to an d rn a are respectively the present age of the universe and the axion mass. For 
the WIMPs, if they decoupled before the temperature was about the electron mass (i.e. if 
the particles (%s) were not affected by the reheating of the photons created by the e + e~ 
annihilations), we find 



where Tr> and m x are respectively the decoupling temperature and the mass of the WIMPs, 
or, if they decoupled after the temperature was about the electron mass (i.e. if the x s were 
affected by the reheating), we find 



Hence, as in the case of light, caustics are expected to be common in the distribution of 
CDM in space. These estimates of primordial velocity dispersions are approximate, but 
in the context of this paper and of galaxy formation in general, 5v a and 5v x are entirely 
negligible. Therefore, the CDM density must be very large at a caustic. Henceforth, unless 
otherwise stated, the velocity dispersion of CDM flow is set equal to zero. 

Detection of the caustic structure in the halos would revolutionize our understanding of 
galaxy formation, and galactic structure. Caustics precisely predict the distribution of dark 
matter in the caustic neighborhood, and hence their gravitational lensing effects [8-11]. Near 
the sharp edges of caustics, we may expect to have special signatures in image magnification 
and in image structure. Point-like background sources and extended sources, such as radio 
jets might thus become probes of elusive CDM particles. Detection of the caustic ring 
structure by gravitational lensing would mean indirect detection the CDM in galactic halos. 
It may even be possible to determine the velocity dispersion of dark matter from lensing 
observations [9], and hence distinguish between axions and neutralinos. Because WIMP 
annihilation rate per unit volume is proportional to the WIMP number density squared, the 
presence of caustics increase the total annihilation rate and change the spatial distribution 
of emission [12] in galactic halos. If photons from WIMP annihilation are observed, the 
caustic locations may be unveiled as sharp lines and hot spots on the sky. 
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In Ref. [9], we studied the density profiles and gravitational lensing effects of outer and 
inner (ring) caustics, in the limit of zero velocity dispersion. For the caustic rings, the 
curvature radii and density depend on the location on the surface of the caustic ring. These 
quantities were expressed as functions of position on the tricusp and estimated using the self 
similar infall model. We presented a formalism which simplified the relevant gravitational 
lensing calculations and applied our formalism to the dark matter caustics in four specific 
cases. In the first three, the line of sight was tangent to a surface where a simple fold 
catastrophe was located. The three cases were distinguished by the curvature of the caustic 
surface at the tangent point in the direction of the line of sight: (i) the surface curves toward 
the side with two extra flows; (ii) the surface curves away from the side with two extra flows; 
and (iii) the surface has zero curvature (a large fraction of the surface of a caustic ring is 
saddle-shaped and therefore has tangent directions along which the curvature of the surface 
vanishes). In the fourth case (iv) studied, the line of sight was at a specific location close to 
a dual cusp and parallel to the galactic plane of the caustic ring. We estimated the lensing 
effects in all these four cases at special points on the caustics. For case (i), which we called 
"concave," only the outer caustics were considered. Cases (ii)-(iv) are valid only for the 
caustic rings. We found that each case has characteristic lensing signatures. In three of 
the cases (ii-iv) there are multiple images and infinite magnification of point sources when 
their images merge (in case (i), there are no multiple images). The effects estimated were 
all small. It was found that the typical magnifications and image distortions are of order 
one % to a few %. The most promising possibility seemed monitoring the sky for distortions 
caused by the dark matter caustics in the images of the extended sources, such as radio jets. 

In this paper, we point out that the lensing effects are enhanced significantly near the 
cusps, in particular, near the outer cusp. We first reparameterize the flow equations near a 
caustic ring. The reparametrization is useful in studying the characteristics, stability and 
gravitational lensing effects of the caustic rings. We express the caustic rings generated by 
the CDM flow as single valued functions z{jp) and p(ip), where z and p are transverse space 
coordinates and ip is an angular variable, parameterizing the tricusp cross-section. These 
equations show that the reason for the existence of the outer cusp which lies in the galactic 
plane, is the same as the other two cusps: both and J| vanish simultaneously at the 
cusps. The parametrization given in [6] is inefficient to reveal this fact, because, the plot of 
the caustic equations do not traverse through the outer cusp. They plot half of the tricusp; 
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the other half is obtained by the reflection symmetry, due to the double valuedness of z(ip), 
in that parametrization. Generically, caustics are surfaces separating regions with differing 
number of flows. One side of a caustic surface has two more flows than the other. The 
density profile near a caustic surface is d = -4^©(<t), where A is called the fold coefficient, 
and a is the distance to the surface, on the side with two extra flows. In the limit of zero 
velocity dispersion, the density diverges as Near the surface of a caustic ring we express 
the fold coefficient as a function of ip. We find that A(ip) is minimum at the middle locations 
(where if) — |,7r, -y) between the cusps. It increases monotonically as one approaches to 
cusps where it diverges. We also study the differential geometry of the caustic ring surface, 
using the Cartan's Structure equations. We obtain Gaussian, mean and principal curvatures 
as functions of ip at any location other than the cusps. As is shown in the Appendix, the 
principal curvatures are the inverses of the principal radii of the surface, and are needed 
for gravitational lensing applications. Using the equations for the fold coefficient and the 
curvature radius that are derived here, in the expressions for the magnification and image 
shift given in Ref. [9], we obtain the lensing formulae for the line of sights parallel to 
the galactic plane of the caustic ring and near tangent to the surface at a given ip. The 
magnification monotonically increases as the line of sight approaches to the cusps where 
it diverges in the limit of zero velocity dispersion. We consider the concave folds of the 
caustic rings to estimate the magnification and the shift numerically, for the first time in 
this paper. We find that, for the axion caustic rings at cosmo logical distances, when the 
observer's line of sight is tangent to the surface at a sample point near the outer cusps, the 
magnification is about 37%. This effect is about 250 times greater than the effect obtained 
for the concave folds considering the outer caustics, and about 37 times greater than the 
largest effect predicted in Ref. [9] considering the convex folds of the caustic rings. We 
also estimate the magnifications at sample locations near the non-planer cusps for cases (i), 
(ii), (iii), and constrain the expected magnifications at the cusps in the presence of finite 
velocity dispersion. When the velocity dispersion is small but nonzero, the divergence at 
the cusps is cut off, because the caustic surface gets smeared over some distance in space 
and the cusps are smoothed out. The lensing effects are no longer infinite at the cusps, but 
merely large. The lower and upper bounds of the effective velocity dispersions of the axion 
and WIMP flows in galactic halos may be used to constrain the lensing effects at the cusps. 
For a cosmological axion caustic ring, we find that the magnification may range between 
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3% and 2800% at the outer cusp, and between 2% and 46% at the non-planer cusps. For a 
cosmological WIMP caustic ring, on the other hand, we constrain the magnification between 
3% and 28% at the outer cusp, and between 2% and 5% at the non-planer cusps. Because 
the upper bounds for the magnification at the cusps are obtained considering the minimum 
primordial value of the velocity dispersions for the CDM candidates in space where no small 
scale structure has formed, and also because the observer's line of sight may not exactly be 
parallel to the galactic plane of the caustic ring in general, they should be regarded with 
precaution. Moreover, a triangular feature in the Infrared Astronomy Satellite (IRAS) map 
of the Milky Way plane was interpreted [13] as the imprint on baryonic matter of the caustic 
ring of dark matter nearest to us. The nearby caustic ring, the fifth ring of the Milky Way, 
is 1 kpc away (in the direction of observation) from us. We also estimate the magnifications 
for the line of sights close to the cusps of this ring. Unfortunately, even near the outer cusp, 
the gravitational lensing due to a caustic ring only a kpc away from us, is too weak to be 
observed with present instruments. Finally, we study the caustic rings in the Catastrophe 
Theory point of view. We derive the Catastrophe Function of the triaxial caustic rings 
and obtain the flow equations as the equilibrium points of this Catastrophe Function. The 
analysis of the Stability (Hessian) Matrix show that the caustic rings are structurally stable. 

To see how caustics of CDM form in galactic halos, it is instructive to discuss the phase 
space structure of CDM halos. Since the particles are in three-dimensional (3D) space, the 
phase space is six- dimensional. The small velocity dispersion Sv means that CDM particles 
lie on a very thin 3D sheet in 6D phase space. The thickness of the sheet is Sv of the particles. 
The number of CDM particles is enormous in terms of astronomical length scales. Hence, 
the sheet on which the particles lie in phase space is continuous. The flows can be described 
in terms of the evolution of this sheet. When a large overdensity enters the nonlinear regime, 
the particles in the vicinity of the overdensity fall back onto it. This implies [6] that the 
phase space sheet "winds up" in clockwise fashion wherever an overdensity grows in the 
nonlinear regime. Before density perturbations enter the nonlinear regime, there is only one 
value of velocity (i.e., a single flow) at a typical location in physical space, because, the 
phase space sheet covers physical space only once. On the other hand, inside an overdensity 
in the nonlinear regime, the sheet wraps up in phase space, turning clockwise in any two 
dimensional cut (r, r) - - where r denotes the radial coordinate and f is the associated 
velocity — of that space. The outcome of this process is an odd number of flows at any 
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point in a galactic halo. One flow is associated with particles falling through the galaxy for 
the first time (n = 1); another is associated with particles falling through the galaxy for 
the second time (n = 2), and so on. At the boundary between two regions, one of which 
has n flows and the other n + 2 flows, the physical space density is very large, because, the 
phase space sheet has a fold there. At the fold, the phase space sheet is tangent to velocity 
space, and hence, in the limit of zero velocity dispersion (5v = 0), the physical space density 
diverges, since it is the integral of the phase space density over velocity space. The structure 
associated with such a phase space fold is called a "caustic." Because the caustic surfaces 
occur wherever the number of flows changes, they are topologically stable. As mentioned 
earlier, in the limit of zero velocity dispersion, the density diverges as d ~ when the 
caustic is approached from the side with n + 2 flows, where a is the distance to the caustic. 
If the velocity dispersion is nonzero, this divergence is cut off, because, the location of the 
caustic surface gets smeared over some distance 5x. For the dark matter caustics in galactic 
halos, Sx and 8v are related [5] by 

ox ~ , (4) 

v 

where v is the order of magnitude of the velocity of the particles in the flow and R is the 
distance scale over which that flow turns around (i.e., changes its direction). For a galaxy 
like our own, v = 500 km/s and R = 200 kpc are typical orders of magnitude. Using the 
estimates of the primordial velocity dispersions of the CDM candidates (Eqs. 1-3), one finds 
[9, 10] that axion caustics in galactic halos are typically smeared over 

/lO- 5 eV\ 4 ' 7/5 ' 7 

(5x Q ~6-10 4 km , (5) 

\ m a ) 

as a result of their primordial velocity dispersion; whereas WIMP caustics are smeared over 

5:r x ~3-10 10 km(^y /2 , (6) 
if the particles were not affected by the reheating, or 

5x x ~5.10 10 km^y /2 , (7) 

if the particles were affected by the reheating. It should be kept in mind, however, that a 
CDM flow may have an effective velocity dispersion that is larger than its primordial velocity 
dispersion. Effective velocity dispersion occurs when the sheet on which the dark matter 
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particles lie in phase space is wrapped up on scales that are small compared to the galaxy 
as a whole. It is associated with the dumpiness of the dark matter falling onto the galaxy. 
The effective velocity dispersion of a flow may vary from point to point, taking larger values 
where more small scale structure has formed; and taking the minimum primordial value 
where no small scale structure has formed. For a coarse-grained observer, the dark matter 
caustic is smeared over 5x given by Eq. 4 where 5v is the effective velocity dispersion of the 
flow. Little is known about the size of the effective velocity dispersion of dark matter flows 
in galactic halos. However, the sharpness of the triangular feature's edges (in the IRAS map 
of the Milky Way plane) implies [13] an upper limit of 15-20 pc on the distance Sx over 
which that caustic is smeared (and hence an upper limit of order 50 m/s on the effective 
velocity dispersion of the corresponding flow). We will use this upper limit and the lower 
limits Eqs. 5 and 7 for the smearing out of caustics to constrain the density and lensing 
estimates at the cusps of the caustic rings in Sects. Ill and V. 

Evidences for the caustic rings were found [13, 14] in the distribution of bumps in the 
rotation curves of spiral galaxies, including the Milky Way. The self-similar infall model of 
galactic halo formation predicts that the caustic ring radii [5] 

{a „ : „ = 1,2, . . .} * (39, 19,5, 13, 10, 8, , , .)kpc . (£=) (H) , (8) 

where j max is a parameter, with a specific value for each halo, which is proportional to the 
amount of angular momentum that the dark matter particles have [15, 16], h is the present 
Hubble constant in units of 100 km/(s Mpc), and i> ro t is the rotation velocity of the galaxy. 
The infall is called self-similar, if it is time independent, after all distances are rescaled by a 
time-dependent scale R(t) at time t, and all masses are rescaled by the mass M{t) interior to 
R(t). In the case of zero angular momentum and spherical symmetry, the infall is self-similar 
if the initial overdensity profile has the form ^ = (If 1 ) , where M and e are parameters 
[17]. In CDM theories of large scale structure formation, e is expected to be in the range 
0.2 to 0.35 [15, 16]. In that range, the galactic rotation curves predicted by the self-similar 
infall model are flat [17]. Because the caustic rings lie close to the galactic plane, they cause 
bumps in the rotation curve, at the locations of the rings. In a study [14] of 32 extended 
and well-measured external galactic rotation curves, evidence was found for the law given 
in Eq. 8. In the case of the Milky Way [13], the locations of eight sharp rises in the rotation 
curve fit the prediction of the self-similar model at the 3% level. Thus, on theoretical and 
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observational grounds, the caustic ring radii a n {n = 1,2,3..) obey the approximate law 
a n ~ \ jn. Evidence for the first outer caustic (the caustic of second turnaround) has been 
recently found [18] in the NGC 5846 group of galaxies. The plot of the surface number 
density of galaxies was inferred to mark the first outer caustic of radius R\ by an abrupt 
density drop (to a roughly constant value) and a transition from a large velocity dispersion 
interior to Ri, to a very small velocity dispersion exterior to R 1 , as predicted by the CDM 
infall models [17]. 

The outline of this paper is as follows. In Section II, we reparameterize the flow equations 
that were introduced in Ref. [6]. Using this reparametrization, we obtain the caustic ring 
equations associated with the CDM flow, as single valued functions of an angular parameter 
ifj. In Sees. IIIA-B, following the procedure given in Ref. [9], we derive the density on 
the caustic ring. Section IV is devoted to the differential geometry of the caustic ring 
surface. In Sect. V, using the equations we obtained in Ref. [9] and in Sect. IV of this 
paper, we derive gravitational lensing formulae for the line of sights tangent to the caustic 
surface at any point on the surface. We estimate the magnification near the cusps where 
the effect is largest, for the caustic rings at cosmological distances (and also for the nearby 
fifth ring of our own galaxy), in the limit of zero velocity dispersion. We calculate the 
upper and lower bounds for the lensing effects at the cusps in the presence of finite velocity 
dispersion. The reparameterized flow equations are also useful in studying the CDM caustics 
in the Catastrophe Theory point of view. In Sec. VI, we present the correspondence of our 
formulation with the Catastrophe Theory. We derive the Catastrophe Function of the triaxial 
caustic rings and obtain the flow equations as the equilibrium points of this Catastrophe 
Function. The analysis of the Stability (Hessian) Matrix show that the caustic rings are 
structurally stable. In Sec. VII, we summarize our conclusions. 

II. FLOW EQUATIONS AT A CAUSTIC RING 

The caustic rings are closed tubes whose cross-section is a LL 4 catastrophe [5, 6]. They 
are located near where the particles with the most angular momentum in a given inflow are 
at their distance of closest approach to the galactic center. For simplicity, we study caustic 
rings which are axially symmetric about the z-direction, as well as reflection symmetric 
with respect to the z = plane. The dark matter flow is then effectively 2-dimensional. 
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Although the flow is in 3D space, in general, the dimension along the tube is irrelevant as far 
as the caustic properties are concerned. Particles actually move in this direction, however, 
the motion is a simple rotation along the azimuthal direction. To discover the essential 
properties of the caustic tubes, it is sufficient to consider the cross-section perpendicular to 
the trivial direction. Thus, the flow is in effect two dimensional, even in the absence of axial 
symmetry. The essential properties of caustics are invariant under continuous deformations. 
To proceed further, we assume axial symmetry of the flow in this paper. We throw away the 
irrelevant ip coordinate and choose p and z as the cylindrical coordinates of the particles. In 
galactocentric cylindrical coordinates, the flow at such a caustic ring is described [6] by: 

p = a + ^(r-r ) 2 - S -a 2 (9) 
z = bra , (10) 

to lowest order in an expansion in powers of r and a; see Ref. [10] for a detailed derivation. 
The parameters r and a label the particles in the flow, p is distance to the z-axis. r is the 
time a particle crosses the z = plane, a is the declination of the particle, relative to the 
z = plane, when it was at last turnaround. The constants a, b, s, u and t are characteristic 
of the caustic ring. 

In this paper, we first make the following reparametrization: 

Xi = ^ (to - t) , X2 = ^|a, (11) 

so that 

P = a + X i 2 -X2 2 (12) 

z = 2((^- X i)X2, (13) 

where ( = These reparameterized flow equations allow us to obtain caustic equations in 
terms of single valued functions p(jp) and z(tp) which are useful in studying the Differential 
Geometry (Sect. IV), Gravitational Lensing (Sect. V) and Catastrophe Theory (Sect. VI) 
of the caustic rings. 

Let us obtain the number density d(p, z, t) of the CDM flow, near a caustic ring at time 
t. The total number of particles 

N = J pdpdz2%d(p,z,t) , (14) 
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can also be written in parameter space as 

d 3 N 



/d 1\ f d 1\ 

dx\dxid<$> , (xi,X2,(f>) = \dx\dxi , , {XuX-i) ■ (15) 

dxidx 2 d(j) J dxidx2 

In the last equation, changing variables (xi,X2) ~^ (,Pi z )i an d using the transformation 
identity 

d 2 N 



A2j\j n{p, z, t) 



dxidx: 



ixuX2) 



det 



d(xi,X2) 
, d(p,z) 



, (16) 



(Xl, X2)=(X1, X2)j(p, z, t) 



we obtain 



f n(p,z,t) ! 

N — dpdz a a (Xi(p,z,t), X2{p,z,t)) jT — 

J ^ dxidx2 \D 2 {Xi, 



X2) 



(17) 



(Xl, X2)=(X1. X2)j(p, z, t) 



Here, we used the facts that the determinant of the Jacobian of the transformation (xi, X2) 
(p, z) is equal to the reciprocal of the Jacobian of the inverse transformation 

d(p,z) \ 



D 2 (xi,X2) = det 



(18) 



d(Xi,X2)J 

and (xi, X2)j, with j = I. . .n, are the solutions of p = X2; t) and z = zfoi, X2; 0- The 
number of distinct solutions (flows) n is a function of the location in space p, z, and time t. 
Comparing Eqs. 14 and 17, we find the number density of particles in physical space as 



d(p,z,t) 



n(p, z, t) 



d 2 N 



— y 

2vrp p[ dxidx2 



((Xi,X2) j(p, z, t)) 



D 2 (xuX2) 

where the determinant D 2 (xi,X2) of the Jacobian matrix 



, (19) 



(Xl:X2) = (Xl:X2)j(p, Z, t) 



V{Xl,X2) = 


( 9p_ 

dxi 


9p_ \ 

9X2 


= 2 


( X, 


-* \ 




dz_ 

\ 9 X i 


dz 

9X2 ) 




V -Cx2 


C(Vp-xi) J 



(20) 



is 



D 2 (xi,X2) = -4C 



v , v 2 P 



(21) 



Caustics occur wherever D 2 vanishes, i.e. wherever the map (xi,X 2 ) —> (p, z ) is singular. 
Hence, at the caustic, in the limit of zero velocity dispersion (5v = 0), the density diverges. 
The divergence is cut off, if the Sv is nonzero, because, the location of the caustic surface 
gets smeared over some distance Sx in space (Eqs. 5-7). 
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Now, to obtain the single valued caustic equations, notice that, in Eq. 21, the curves for 
which D 2 = are circles (for p^O) centered at (xi, X2) = (-^y, 0) with radius (Here we 
assume that xi an d X2 are Cartesian coordinates, if not, the critical curves are ellipses for 
p 7^ 0). The critical circles degenerate into the singular point (xi,X2) = (0,0) in the limit 
p = 0. The parameter representations of the critical circles can be given as: 

Xi = ^(l±cosV>) » X2 = : YsinV, (22) 

where the angular variable tp G [0, 2n]. Since the CDM density diverges where D 2 vanishes, 
by substituting Eq. 22 into Eqs. 12-13, we find parametric equations describing the cross- 
section of the caustic surface in the (p, ^)-plane: 

p(^) = a + ^cos^(cos^±l) (23) 
p 

Z W) — (7; sin'0(l =F cos^O . (24) 

The parameter a is the radius of the caustic ring and p is the longitudinal dimension of the 
caustic cross-section. The parameter (, as is shown in the below, is proportional to the ratio 
of the latitudinal dimension q to the longitudinal dimension p: ( = Figure 1 shows a 

plot of the cross-section, which we call tricusp. 
The tangent vector field on the cross-section 



— (p(tp)x + z(ip)z) = -- 



( sin 2ip ± sin ipj x — £ ( cos ip =F cos 2 



(25) 



vanishes at ip — 0, -y, -y (^ = y , 7r, y) for the upper (lower) sign choice of Xi- At these 
values, where both and simultaneously vanish, caustic has dual cusps; see also Sect. 
VI. Thus, the reason for the existence of all three of the cusps is the same: simultaneous 
vanishing of and ||. This point can not be seen by the parametrization of Ref. [6] 
because, the tricusp obtained there is made up of two separate pieces on top of each other. 
As the parameter r ranges in the interval [0, r ] the functions p(r) and z{r) of [6] describing 
the caustic cross-section plots only the half of the tricusp. The other half is obtained by 
the reflection symmetry of the flow equations, or, i.e. by the double valuedness of the 
caustic equation z(t). Neither p(r), nor z(t) traverses through the outer cusp. By the 
new parametrization, on the other hand, as the angular parameter ip ranges in the interval 
[0,27r], the functions p(ip) and z(ip) (Eqs. 23-24) plot the whole tricusp as a single piece, 
traversing through all the cusps, including the outer one. 
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Notice also that, the vanishings of both the determinant D 2 (Eq. 21) and tangent vector 
t (Eq. 25) at three specific values of ip, are independent of (. This means that, the existence 
of the caustic rings is independent of the ratio of transverse dimensions p and q (i.e. the 
parameter Q, and the number of cusps is always three. For £ = 1, the map ip — > ip + ^ 
transforms the three cusps of the caustic into one another. Thus, in this case, the tricusp 
has a Z 3 symmetry. In the language of Catastrophe Theory, the tri-axial tricusp, where the 
three cusps make up an equilateral triangle (£ = 1), is a -D_ 4 catastrophe. To avoid the 
clutter of the sign choices, we adopt the upper sign conventions in Eqs. 23-24 and use 

p p 
p(ip) = a + -cos + cos ip) , z(ip) — (-sin — cos tp) , (26) 

in the rest of the paper. Since both p(ijj) and z{jp) are extremized for tp — 0, 4^, we 
can deduce the dimensions of the tricusp by evaluating p{jp) and z{ip) at those locations: 
(p(0),*(0)) = (a+p,0), (p( 2 f),z(f)) = (a-f,C^#p),and (p(f )) = (a-f,-^). 
Thus, as is mentioned before, the dimensions of the cross section in the p and z directions 
are p and q = C^4%>, respectively. In terms of p and q, the locations of the three cusps in 
physical (p, z)-plane are therefore: (a +p, 0), (a — |, |), and (a — |, — |). 

In physical space, caustics are located where the number of flows changes by two. Inside 
(outside) the tricusp, there are four (two) flows. To count the distinct flows, let us restrict 
ourselves to the z = plane, where, from Eq. 13, we either have (i) Xi — \fp~-, or > (h) X2 = 0. 
If (i) is the case, then Eq. 12 gives p = p — X2 2 < Po- Therefore, at least two distinct flows 
exist in the region p < p Q . They are parameterized as (xi, X2) = (yfp, i^Po — p)- For both 
of the flows, the velocities in the p direction are the same: 

We use ~ because the flow is not stationary, (if one assumes that the flow is stationary, then 
~ is replaced by =). The velocities of the flows in the z direction 

m~-Tr = • (28) 

are opposite to each other. For the flow (xi, X2) = (y/P, V Po — p), ff — —C\J^ u {po — P) < 0, 
thus we call it the "down-flow." For the flow (xi, X2) = (\/P, ~VPo ~ P), 011 the other hand, 
|| = (^j2u(po — p) > 0, thus we call it the "up-flow." From Eq. 19, we know that the 
contribution of each of the flows to the density is inversely proportional to the absolute 
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value of D 2 = 4£(p — Po) = 4£(p — a ~ p)- Near p = a, as p — a — > ± , -D 2 — > ~^(P- 
Therefore, near p = a, each of the contributions of the down and up flows is finite: d oc 
If, on the other hand, one approaches to the outer cusp from the inside of the tricusp, as 
p — po —* 0-, (recall that, the down and up flows exist only where p < po), the densities of the 
up and down flows diverge: d ~ ^ ( n ) i s the case, then Eq. 12 gives p = a + Xi 2 > a - 

Therefore, in the region p > a, at least two distinct flows exist. They are parameterized as 
(Xi> X2) = (±VP — a i 0)- F° r both of the flows, the velocities in the z direction are zero: 

Y t --% = -^CX 2 = . (29) 
The velocities of the flows in the p direction 



are opposite to each other. For the flow (X15X2) = (Vp~-~a, 0), |f = \J^u(p — a) > 0, 
thus, we call it the "out-flow." For the flow (xi,X2) — (— \J P — a, 0), on the other hand, 
|^ = —^J2u(p — a) < 0, thus, we call it the "in-flow." For in and out flows, D 2 = — 4^(p — 



a =f yp(p — a)), respectively. Near p = a, as p — a — > 0+ (recall that, the out and in flows 
exist only where p > a), D 2 — > ±4^y / p(p — a). Therefore, near p = a, the densities for the 
out and in flows diverge: d ~ ^ — . 

Next, let us analyze the density properties near the dual cusp, where the behavior depends 
upon the direction of approach. For z = and p — p — > 0_ (inside the caustic tube), 
there are four flows: in, out, up and down. For p — po — > + (outside the caustic tube), 
there are two flows: in and out. Here, let us use this opportunity to point out that the 
density of the out-flow diverges in both sides of the cusp, whereas that of the in-flow remains 
finite. To see this, let us express the radial coordinate as p = a + p + e, where e — > 0±, 
depending on the side we approach to the cusp. Therefore, we can express the determinant 
as D 2 ~ — 4C[p + e =F p(l + j-)], for the out and in flows, respectively. Thus, for the out 
flow, D 2 ~ — 2(e, and hence the density diverges, d ~ A, in both sides of the cusp. For 
the in flow D 2 ~ — 8(p, thus, the density of the in flow remains finite in the both sides. 
Notice that, near the cusp, unlike the region p — a < where both of the down and up 
flows remain finite as one approaches the fold caustic from the outside, there is a flow (the 
out flow) whose density diverges as one approaches the caustic tube from the outside of the 
cusp. Hence, the cusps are more divergent then the folds. 
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Finally, for the sake of completeness, let us calculate the density near the cusp, again at 
p = p , but now, as z — > 0±, in terms of our new parameters. Since p = a + Xi 2 — X2 2 , at 
p = Po = a + p, we have 

X2 2 = (xi + Vp)(Xi~Vp)- (31) 
Using Eqs. 31 and 13, the determinant D 2 (Eq. 21) can be written as 

D 2 = -4C [2(xi - VP) 2 + 3VP(Xi - VP)} = -7 f — ) + • (32) 

We also infer from Eqs. 13 and 31 that, near the galactic plane, at p = p , z — > — — > 0, 
because the parameters (xi,%2) - ^ (>/P> 0) there. Therefore, in this limit, the second order 
term ( — ) 2 in Eq. 32 can be neglected next to the first order term — . Hence, we have 
D 2 — —6 [Cpz 2 ] 1 ^ 3 for one of the flows, and the density d ~ -^73 f° r this flow. Thus, the 
density also diverges in this particular limit near the cusp. 

Note also that, the tube caustic collapses to a line caustic in the limit p — > with ( fixed. 



In this limit, Eq. 21 gives D 2 = -4CVXi 4 + X2 4 + 2xi 2 X2 2 = ~^C 2 (p ~ a) 2 + z 2 ■ Hence, 
the density 

1 d 2 N 1 

d(p, Z) = : : ; , (33) 

8np d X ld X 2 ^2( p _ fl )2 + Z 2 

in the limit p — > 0. 



III. DENSITY AT A CAUSTIC RING 



In the previous section, we have restricted ourselves to the p-axis and worked out the 
density of the flows near the points p = a and p = p Q . They are the points where the 
parameters that describe the flow, cross the degenerate critical circle xl + X 2 ~ \/VX\ — 
on the p-axis (see also Sect. VI). These locations are interesting, because, as we will see in 
Sect. VI, on the p-axis at p = a a simple fold (A 2 ) catastrophe, and at p = p a dual cusp 
(A_ 3 ) catastrophe, necessarily occur. In this section, using the equations that describe the 
caustic ring in terms of the angular parameter vp, and following the procedure given in Ref. 
[9], we derive the density profile near the caustic surface. In Sect. Ill A, we derive a density 
formula valid everywhere except at the cusps (our mathematical formulation is invalid at 
the cusps). In the limit of zero velocity dispersion (5v = 0), the density diverges when one 
approaches a caustic surface, on the side which has two extra flows, as the inverse square 
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root of the distance to the surface. If the velocity dispersion is small, but nonzero, this 
divergence is cut off, because, the location of the caustic gets smeared out. So the density 
at the caustic is no longer infinite, but merely very large. In Sect. Ill B, we deal with the 
density near the cusps. Our results are applied in Sect. V, to derive the gravitational lensing 
signatures of the caustic rings. 



A. Density everywhere except the cusps 



We choose an arbitrary point on the surface of the caustic ring parameterized by 
(Xi(^*)i X2{i>*)) where XiW an d X2C0) are given by Eq. 22. We assume that the point is 
not at one of the three dual cusps located at ip — 0, and 4?, and use ip* to indicate this 



exclusion. Whenever ip, without a subscript "*" is used, the dual cusps are included. The 
physical space coordinates (p*, z*) = (p(i/)*), 2(1/)*)) are given in terms of ■?/>* by Eqs. 26. The 
vanishing of Z^^*) = detV (■?/>*) implies the existence of a zero eigenvector of the matrix 
X>('0*) and a zero eigenvector of its transpose 

/ 

V T (^) = 2 



(34) 



Let us define 9* = 9 (■0*) such that 



sinf?* 

cos 9* 



. 



Hence, we have 

Xi(^*) sin #* - CX2(^*) c °s 9* = 
~X2(ip*) sin^ + C(VP - Xi(fp*)) cos0 * = • 
Equations 36 and 37 imply 

±ClX2(^)l 



sin 9* 



cos 9* 



±((Vp-xM) 



_ ±C 

v/xi(^) 2 + c 2 x2(^) 2 ^X2(^y + e(vp-xi(^)) 2 ^Je + cot 2 f 

±Xl(V>*) ±1X2(^)1 ±1 



C 2 tan 2 & 



(35) 

(36) 
(37) 

(38) 
(39) 



/xi(^) 2 + C 2 x 2 (^) 2 v / x 2 (^) 2 + C 2 ( v / P-xi(^)) 2 V 1 

Now, let us define a new pair of Cartesian coordinates (r/, a) related to (p — p*, z — 2*) by a 
rotation of angle 0* + |: 

' ~- — W 

(40) 




cos 9* — sin 
— sin 6L — cos 9* 



P~ P* 
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We now show that a is the coordinate in the direction orthogonal to the caustic surface at 
(p*, z*). Consider small deviations about (xi^*)? X2(i J *)) m parameter space: (xi,X2) — 
(Xify*) + Axi,X2{ip*) + A% 2 )- Equations 12, 13 and 20 imply 

/ 



\ Az ; 




0(A X i 2 ,Ax2 2 ,AxiAx 2 ) , 



(41) 



where 



Ap = P ~ P* = 2(xi(V>*) Axi - X2M Ax 2 ) + 0(A X i 2 , A X2 \ A X iAx 2 ) 

Az = z - z* = 2C(-x 2 (^)A X i + (Vp~ Xi(^*))Ax 2 ) + 0(A X i 2 , A X2 2 , A Xl A X2 ) . (42) 



The expansion of a, defined by Eq. 40, in powers of Axi and A% 2 yields 



( 7 = 0(Axi 2 ,Ax 2 2 ,AxiA X2 ) , 



(43) 



because, due to Eq. 35, the first order terms vanish. The fact that a is second order in Axi 
and Ax2 shows that a is the coordinate in the direction perpendicular to the caustic surface. 
This can be seen as follows. The scalar product of the tangent vector dt = pdp + zdz and 
the normal vector h = pn p + zn z at ip* is zero. Thus, at using Eq. 41, in first order, we 
have 

/ 



dt-n=(Ap Az) 



n. 



n. 



(A X i A X2 )D 1 (t 



n r 



0, 



(44) 



n. 



which implies D T {il)^)n = 0. Recall, however, that we defined the zero eigenvector, and 
hence 9(ip*), such that Eq. 35 is satisfied. Therefore, the zero eigenvector must be the 
normal vector n = psin#* + £cos#*, where 9* is the angle between the p-axis and the 
tangent vector at (p*,z*); see Fig. 2. Thus, in first order, a is parallel to the normal. We 
choose its direction pointing inward, therefore, at the given of Fig. 2 we have 



a = —p sin 9* — z cos 9* . 



(45) 



To obtain the density profile of the caustic ring near the point under consideration, we need 
D 2 (ri,a) to order y/a. So, we calculate a to second order in powers of Axi and Ax 2 , and 
D 2 and 77 to first order. Using Eq. 40 in the expansion: 

1 



a = 



A X i 2 + 2— — — 

i>* OX1OX2 



dxi 2 



A A d2 <? 

AxiA X2 + j— 2 



A X 2 2 



(46) 
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FIG. 2: An arbitrary point on the tricusp is labeled by ip*. Its physical coordinates are (p*, z*). A 
new Cartesian coordinate system (a, rj) is defined there such that a is perpendicular to the caustic 
surface. It is rotated relative to the (p, z) coordinates by an angle 0(Y>*) + |. 



we find 



o = (Axi Ax 2 ) 



' - sin 6>* C cos 0* \ ( Axi 
^ C cos sin ^* y y 



(47) 



Equation 21 yields 
£ 2 



dxi 
= -8C 



V>* dx2 



AX2 + 0(Axi 2 ,Ax2 2 ,AxiAx 2 ) 



Xi(^*) - \ J Axi + X2(^*)Ax 2 



0(A Xl 2 ,A X2 2 ,A Xl A X2 ). (48) 



To find 7], we use Eqs. 42 in Eq. 40: 



drj 



, A * 1 + Tr- 
ip* C)X2 



A X2 + 0(A X i 2 , A X2 2 , AxiA X2 ) 



= 2 [xi {ip* ) cos 0* + CX2 (^* ) sin 0,] Axi - 2 [x 2 (^* ) c °s + C ( y/p - Xi O* ) ) sin 0*] A X 2 ■ (49) 

Equation 49 can be recast using Eqs. 38 and 39 as 

V = 2 V / Xi 2 (^) + C 2 X2 2 (^)Axi - 2 V / X2 2 (^)+C 2 ( v / P-Xi(^)) 2 Ax2 • (50) 

The fact that D 2 (ip*) = implies ^/p — Xi(V ; *) — X2 2 (V ; *)/Xi(V ; *)- We can, therefore, write 
Eqs. 48 and 50, at any ip*, as a matrix equation: 

\ 

(51) 



D2 

8C 



\ 



V 2 v /xi 2 W*)K 2 X2 2 W*) / 



V xi(^*) — |X2(-0*)I / 



Axi 

V A X2 y 
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Equation 51 can be inverted to obtain Axi and A%2 as functions of D 2 and rj. When the 
result is inserted into Eq. 47, we obtain 



Cxi(^*) (ft 



v / PXiP/'*)'? 2 



4C 2 xi 2 (V>*)+C 2 X2 2 (V>*) 



16 (xi(^) - f ) lx 2 (^)l\/xi 2 (^) + C 2 X2 2 (^) ' 



(52) 



This gives 



|Afa,<7)|: 



A 



8C 



pA/(l + cos^) 2 +C 2 sin ip* 



'1+2 cos^*)tan 



a + 



(V 2 



l + C 2 + (l-C 2 )cos^ 



(53) 



Along the a direction (rj — 0), we have 



|£> 2 (<r)| = 2^CP^C(^), 



where 



C(<M = 



\ 



'1 + 2 cos ■0*) tan 



0* 



1 + cos-0*) 2 + C 2 sin 2 i J * ■ 



(54) 



(55) 



The derivative of C(ijj*) implies that, near the caustic surface, the density decreases as one 
approaches to the central points between the cusps where ■?/>* = |,7r, and whereas it 
increases as one approaches to the cusps where ip = 0, and -y. Combining Eqs. 19 and 
54, and minding the factor of two because two flows contribute, we obtain the density profile 
near the surface along the a direction, as 



with 



where we used 



V " 

d 2 M /2C cosa(^ 



d 2 M 



d 2 M 2 d 2 M2( d 2 M4ir( 



cos a (■?/>*) • 



(56) 



(57) 



(58) 



dXidX2 dadr y/us dadr b dfldt b 

Here, we define ^| = dadtl!fcos a as ^ ne mass falling in, per unit time and unit solid angle. 
Note that, A(ip) diverges at each of the three dual cusps because C vanishes there. Moreover, 
one has 



d 2 M 
dVtdt 



= fnV r 



y rot 



4ttG ' 



(59) 
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where v n is the velocity of the particles in the nth caustic ring, and the dimensionless 
coefficients f n characterize the density of the nth in and out flows. In the self-similar model 
[5], we have 

{/„ : n = 1, 2, . . .} ~ (13, 5.5, 3.5, 2.5, 2, . . .) • 1CT 2 (60) 



for e = 0.2. The 



dM 



does not depend sharply upon e. Our estimates are for e = 0.2 



dQdt 

because they can most readily be obtained from Ref. [5] in that case. Combining Eqs. 57 
and 59, we find 



AM = ^% - ^ J— cosa^O-W.) , (61) 
4ir& a n b n V Pn 



where 



[1 + ^COS^*(l + COS l\)-if 



j: UU \ — — : ££n : 1 __J (Q2) 

J\ (1 + 2 cos tan f | [(1 + cos ^) 2 + Q sin 2 ^} I 

Let us note, as a check, that C n (n) = y/2£ n , and hence A n (n) reduce to A 0:n of Ref. [9], which 

are the fold coefficients evaluated at (p, z) = (a n , 0). (The A n are minimum at ^ = f , tt, ^f, 

they increase as one approaches to the cusps where V* = 0,^,^). It was shown in Ref. [6] 

that b n and v n are of the same order of magnitude. Moreover, the ten rises in the rotation 

curve of the Milky Way were interpreted [13] as the effect of caustic rings. In that case, the 

widths p n of caustic rings are determined from the observed widths of the rises. Typically 

one finds p n ~ 0.1 a n . Using this, v n ~ b n and cos «(■?/>*) ~ 1 in Eq. 61, yields 

{4.(*):»=1,2, . .}~(5,6,6,7,8„ . .)10-5^f(^) \±) ''(^t) ■ (63) 

Cm 2 kpC2 V^max/ V - 7 / V 220km / S / 

Let us now estimate the fold coefficient more precisely for the fifth caustic ring of the Milky 
Way, which is the one closest to us and believed to be the most constrained by observation. 
For that ring a 5 = 8.31 kpc, b 5 = 516km/s (if r < 0) or 657km/s (if r > 0), p 5 = 0.134 
kpc, q 5 = 0.2 kpc, f 5 = 0.02, v 5 = 480 km/s and ( 5 = = 1.15. Therefore, if r < 



^) = 1 -™- 10 - 3 7^i{^) ' (64) 



or, if r > 



.gJ^^Ogr / v Iot 



AM = 1.38 • 10" , OT (65) 

cm 2 kpc 2 V220km/sy 

where we have taken cos(a(^*)) = 1. If we expand as a power series around the outer 
cusp, where ^ = i? or ^ = 2ir — we find 



, 3 

. 7 , 



iH 4-C 2 + 16-C 2 ^ f 

T (# )=T n (2tt-i? ) = -=7- + - ^— ^ aJ +G>(fl*). (66) 

V ; V ; ^3 1 + ^ 16^1 + p - 2 
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The points ip* = $ and ^* = 27r — are located at 

p„(tf) = Pri (2vr - i?) = a n + p„(l - h 2 + A# + 0( ^ 6)) 

z n (tf) = -z n {2it - 0) = - ^ + 0(tf 7 )) . (67) 

_i _i 
We see that, as ip* — > 0, ^(ip*) diverges as ip* 2 (for p n = 0.1 a n , J r (^) ~ 0.525 $ z, for 

p 5 = 0.016 a 5 , JF" 5 (^) ~ 0.568 However, as is discussed in Sect. I, the divergence is 

cut off, because, the location of the caustic surface gets smeared out over some distance Sx 

(Eqs. 5-7) and the cusps are smoothed out. Recall, however, that these distance scales are 

estimated taking only the primordial velocity dispersions (Eqs. 1-3) into account. Therefore, 

they provide theoretical lower bounds for 5x. CDM flow may have an effective velocity 

dispersion 5v c s which is larger than its primordial velocity dispersion. Unfortunately, very 

little is known about the size of the 5v e g of CDM flows in galactic halos. As is noted in 

Ref. [9], the sharpness of the edges of the triangular feature in the IRAS map [13] implies a 

rough upper bound of 15-20 pc on 5x. We use this upper bound and the lower bounds for 

5xs, which are calculated using the primordial velocity dispersions of axions and WIMPs, 

to estimate the upper and lower limits for the density and magnification at the smoothed 

dual cusps, in the presence of finite velocity dispersion. To find an upper bound, we first 

assume that the velocity dispersion of the particles are zero, and then estimate the effects at 

a location on the surface, whose transverse distances (Ap, Az) to the cusp are greater than 

the distance that a caustic surface would smear out, if the particles had a finite primordial 

velocity dispersion. The lower bounds for the effects are estimated at a location whose 

distance to the cusp is about the maximum smearing out distance of the caustic surface 

implied by the triangular feature in the IRAS map of the Milky Way. Thus, to estimate 

the lower limits of the density and magnification for the smoothed cusps, we choose the 

locations that are A^* = radian away from the locations that the cusps would occur if 

the velocity dispersion were zero. Similarly, the upper bounds for the density and lensing 

effects at the smoothed cusps are estimated at a location on the surface, whose transverse 

distances (Ap, Az) to the cusp are greater than the minimum distance that a caustic surface 

would smear out, if the particles had only the primordial velocity dispersion. We choose the 

locations that are respectively Aip* = and A^* = ^ radian away from the cusps of the 

axion and WIMP caustic rings, to estimate the upper limits. 

Thus, the lower bounds for the density and magnification at the smoothed cusps are 
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obtained at ^* = ±^,^±^,^±7^. The locations of these points with respect to the 
nearby cusps, in the limit Sv = 0, are: (Ap n , Az n )(±j^) ~ (1.26, ±1.76 • lCT 1 ^) • lO^Vn, 
(Ap n , Az n )(f -^5) = (Ap n , Az„)(f ±^) ~ (7.82 • 10" 1 , C„) • lO" 1 ^, and (A Pn , Az n )(f + 
f^) = (Ap n ,Az n )(f - ^) ~ (4.78 • 10- 1 , 1.18C™) • 10" V- We have, for the fifth ring, 
(Ap 5 ,Az 5 )(±^) ~ (16.9pc,±2.7pc), (Ap 5 , Az 5 )(f - ^) = (Ap 5 , Az 5 )(f + ^) ~ 
(10.48 pc,15.44pc), (Ap 5 ,Az 5 )(f + ^) = (Ap 5 , Az 5 )(f - ^) ^ (6.40 pc, 18.15 pc), 
which are all about the size of the upper limit for smearing out of the cusp locations 
estimated using the triangular feature in the IRAS map of the Milky Way. To es- 
timate the upper bounds for the density and magnification near the cusps of the ax- 
ion caustic rings, we choose = ±7^5, ± 7^9, ± j§^- The locations of these 
points, with respect to the nearby cusps, in the limit Sv = 0, are: (Ap n , Az n )(±j^) ~ 
(1.32, ±1.84 • 10- 4 Cn) • 10- 7 p n , (Ap n ,Az n )(f ± ^) ~ (6.58 • 10" 1 , l.UQ ■ 10~ 7 p n . For 
the fifth ring (Ap 5 , Az 5 )(±^) ~ (1.76 • 10- 5 pc,±2.83 • 10- 9 pc), (Ap 5 , Az 5 )(f ± 7^0) = 
(Ap 5 , A^ 5 )(^ =F ^) ~ (8.82 • 10~ 6 pc, 1.75 • 10~ 5 pc). These distances are all grater than 
the lower bound for smearing out Sx a of the caustic surface location, estimated using the 
primordial velocity dispersion of the axions. For the WIMP caustic rings, considering Eq. 
7 (i.e. assuming the particles were affected by the reheating), the upper bounds for the 
density and magnification are estimated at ^ = ±t^,^±^,4^±^. The locations 
of these points with respect to the nearby cusps are: (Ap n , Az„)(±|g) ~ (1.32, ±1.84 • 
KT 2 C„) • 10- 3 p„, (Ap n ,Az n )(f + %) = (A Pn ,Az n )(f -#) ^ (6A2.10-\1.15Cn)-10- 3 Pn, 
(Ap n , Az n )(f - f 5 ) = (Ap n ,Az n )(f ± ~ (6.74 • 10" 1 , 1.13C„) • 10" 3 p„- For the fifth 
ring, we have (Ap 5 , Az 5 )(±f- 5 ) ~ (1.76 • 10- 1 pc,±2.83 • 10- 3 pc), (Ap 5 , Az 5 )(f ± #) = 
(Ap 5 , Az 5 )(f - f 5 ) ~ (8.60 • 10- 2 pc, 1.77 • 10" 1 pc), (Ap 5 , Az 5 )(f - = (Ap 5 , Az 5 )(f + 
jg) ~ (9.03 • 10~ 2 pc, 1.74 • 10 _1 pc). These distances are all grater than the lower bound 
Sx x for smearing out of the caustic surface location, estimated using the primordial velocity 
dispersion of the WIMPs. 

We, therefore, estimate the lower bound of the T at the outer cusps of the caustic rings, 
assuming p n = 0.1 a n , ( n = 1, as jF(±y^) ~ 0.85 (for the fifth ring of the Milky Way, 
jF 5 (±y^) ~ 0.9). To predict the upper bounds for the T near the cusps of the caustic rings, 
we distinguish between the axion and WIMP caustic rings. At the outer cusp of a caustic 
ring, assuming a n = 0.1 a n and ( n = 1, we find an upper bound as jF(^y) ~ 25.6 for the 
axion rings. If the fifth caustic ring of the Milky Way (where p 5 = 0.016 a 5 , ( 5 ~ 1.15) is 
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axionic, then the upper bound is jF 5 (^y) ~ 27.8. Thus, the function T n in Eq. 63 may 
enhance the density at the outer cusp of an axion caustic ring in the range between 0.9 and 
26 times, depending on the size of the effective velocity dispersion. Taking a sample point at 

■0* — ^ — ifo an< ^ assumm g Pn — 0.1 a n , ( n = 1, we find = -^(ifif ) — 3 ' (^ or ^ ne ^h 

ring of the Milky Way ^5(3^) — 3.2). Therefore, at the point tp* = whose location with 
respect to the outer cusp is (Ap n , Az n ) = (7.4- 1CT 4 , 7.75- 10~ 6 ( n )p n , the factor T triples the 
fold coefficient. For the fifth caustic ring of the Milky Way, the location given by ip* = is 
9.92 ■ lCT 2 pc (about 3 ■ 10 12 km) away from the location of the cusp that would occur in the 
limit Sv = 0. The fold coefficient at this location A 5 (j^) ~ 5 ■ 10~ 3 gr/cm 2 kpc^, where we 
averaged Eqs. 64 and 65. For the WIMP caustic rings, the upper bound is JF(±^) ~ 2.6 (if 
the fifth ring of the Milky Way is a WIMP caustic, then the upper bound is JF 5 (±|^) ~ 2.8). 
Therefore, the function T in Eq. 63 may enhance the density at the outer cusp of a WIMP 
caustic ring in the range between 0.9 and 2.6 times, depending on the size of the effective 
velocity dispersion. 

Next, to estimate the fold coefficient at points around the other two non-planer cusps, 
which are located at 

Pn(f ± <>) = Pn(f T$) = a n -^(l-3$ 2 ±V3^+ 3 -^ + 0(#)) (68) 



± <>) = ~z n {% T #) = UvJ^- ( 1 - ^ 2 T i + ? + 0(^)) , (69) 



3 ' " v 3 8 V V27 4 

where we expand T n as: 

3 3 ^(8-^(1 + 3^ (8-^)(l + 3C 2 )i 

For p„ = 0.1 c^, (n = 1, we find F(?f ±&) = J-(^- =F 1?) ^ 0.585 _ 5, which diverges as one 
approaches to the cusps. Near the non-planer cusps of the rings (taking p n — 0.1 a n , ( n — 1) 
we find the upper bound of for the axion caustic rings as F(^±i^) = ^(^zp^j) ~ 28.6 
(if the fifth ring of the Milky Way is an axion caustic, then J-^^f- ± 7500) — -^("f" -F 7500) — 
26.8). Thus, depending on the size of the velocity dispersion, the factor T in Eq. 63 may 
enhance the density in the range between 0.9 and 28.6 times, near the non-planer cusps of 
the axion caustics. For the WIMP caustic rings, on the other hand, the upper bound of T 
near the non-planer cusps, is ± ^) = T(*f =p ^5) — 2.9 (if the fifth ring of the Milky 
Way is a WIMP caustic, then T h i^f ± ^) = =F - 2.7 is the upper limit of T at the 



25 



non-planer cusps). Therefore, the factor T in Eq. 63 may enhance the density in the range 
between 0.9 and 2.9 times, near the non-planer cusps of the WIMP caustics. The locations 
that differ by ±^ from the non-planer cusps of the fifth ring of the Milky Way, are about 
0.18 pc away from the locations where the cusps would occur if the velocity dispersion of 
the flow were zero. 



B. Density at a dual cusp 

In this section, following the procedure given in [9] and assuming that the velocity dis- 
persion 8v = 0, we derive the dark matter density profile near a dual cusp, in terms of the 
new parameters. We choose the outer cusp which lies in the z = plane at p = a + p = p . 
Close to the z = 0-plane, from Eq. 13, we have either xi — \[Vi or X2 — 0, or both. Near 
p = po, however, from Eq. 12, we have xi 2 ~ X2 2 — P- Therefore, near the outer cusp, we 
must have both xi — \/V an d X2 — 0. Defining new dimensionless quantities 

A = — , T=l-2-, X^P^, Z=±, (71) 

y/P y/p P (P 

we recast Eqs. 12-13 as: 

X = -2T + T 2 — A 2 , Z = 2AT, D 2 (A, T) = 4(p (t — T 2 — A 2 ) . (72) 

Close to the cusp, where xi — \fp~i we ma y neglect the terms of order T 2 in Eqs. 72. The 
term of order A 2 cannot be neglected. Therefore, near the cusp, we have 

X = -2T -A 2 , Z = 2AT , D 2 (A, T) = 4(p (t - A 2 ) . (73) 

Inserting T = ^ into the equation for X, one obtains the third order polynomial equation: 

A 3 + XA + Z = . (74) 

The discriminant is: 5 = ) + (y) . If 5 > 0, the cubic equation has one real root, and 
two complex roots which are complex conjugate of each other. If 5 < 0, all the roots are 
real and unequal. For 5 = 0, all the roots are real and at least two are equal. The number 
of real roots is the number of flows at a given location. The tricusp has two flows outside 
and four inside. In the neighborhood of a cusp, however, one of the flows of the tricusp is 
nonsingular and does not participate in the cusp caustic. To include the root corresponding 
to the nonsingular flow near (p, z) = (p , 0), one must keep the terms of order T 2 in Eq. 72. 
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The equation for the caustic surface in physical space is 5 — 0. Indeed, Eqs. 73 imply 
that at the dual cusp caustic, where D 2 vanishes, T = A 2 (therefore, Z = 2A 3 = 2T 3 / 2 and 
X = —3 A 2 = — 3T at the cusp). The discriminant, using Eqs. 73, can also be written as 
5 = 10 ^p (^4 4 + 7A 2 T — 8T 2 ). Thus D 2 = implies 5 = 0, however, the converse is not true: 
5 = does not necessarily imply D 2 = 0, because not all flows at the location of the caustic 
surface are singular. 

Inserting the D 2 obtained in Eq. 73 into Eq. 19, yields 

d _ _1 1 d 2 M " 1 _ 1 1 d 2 M " 1 

where the sum is over the flows (i.e., the real roots of the cubic polynomial Eq. 74). In the 
last equality we used the equation for X (Eq. 73), to rewrite T — A 2 = — \{X + 3A 2 ). If the 
discriminant 5 = (Z/2) 2 + (X/3) 3 > 0, the one real root is 

A=(-f + v / 5) 1 /3_(| + V5)i/3. (76) 

This describes the one flow outside the cusp. Inserting Eq. 76 in Eq. 75, we obtain 

d= J_J_J^_ 1 (7?) 

4ir(p pdxid X 2 \X - 3(-f + V5) 2 / 3 - 3(f + v^) 2 / 3 | ' 

Next, we calculate the density inside the cusp, where 5 < 0. The three real roots of the 
polynomial Eq. 74 are [9] 



A 1 = 2 ] J-^cose, A 2 = 2 ] J-^cos(6 + ^), A 3 = 2^cos(^ + ^) , (78) 

7 / o \ 3/2 

where cos 3$ = — ^ (— ^) and < ^ < |. Inserting them into Eq. 75, and adding the 
individual flow densities, we find 

1 1 d 2 M 1 1 

d = (79) 

2tt( p pdxidx2 \X\ ( v / 3-tan^)sin2^ ' 

Note that using Eq. 58 and cos a ~ 1 reproduces the results obtained in Ref. [9]. 



IV. DIFFERENTIAL GEOMETRY OF CAUSTIC RINGS 



In this section, differential geometry of the axially symmetric caustic rings is studied. 
We exclude the cusps from the discussion and apply the method of moving frames through 
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out the section. We calculate the local quantities, such as metric, Gaussian, mean, and 
the principal curvatures of the caustic surface. The principle curvature radii, which are the 
inverses of the principal curvatures (see the Appendix), are used in the lensing applications 
in Sect. V. 

The caustic surface is completely described by the three vector 

X(ip, 0) = p(ip) cos x + p(ip) sin y + z(ip) z , (80) 

where -0 and are angular variables ranging in the interval [0, 27r], and (x,y,z) are the 
Cartesian unit vectors. The functions p(ip) and z(tp) are given in Eq. 26 which we copy here 

P P 
p(ip) = a H — cos -0(1 + cos-0) , z(ip) = (- sin -0(1 — cos-0) . 

At each point P of the caustic surface, excluding the dual cusps where -0 = 0,^, and 
-y (denoting this exclusion by * as in the previous section), we choose three differentiable 
orthonormal vectors e^, % = -0*, 0, n, with the inner product 

ii ■ ij = 5ij , (81) 

such that e n is the unit normal h, and hence and span the tangent plane: 
A _ X^ t p' (cos (px + sin <py) + z'z 

e *- = pTJ\ = VFT^ (82) 

— * 

= — = — sin x + cos y (83) 
ll^ll 

- z' (cos (j)x + sin (f)y) + p'z 
e n = n = e^xe 4> = = ^ . (84) 

The subscript comma "," indicates partial differentiation with respect to the index that 
follows it. The prime denotes derivative with respect to -0*. The reason for the exclusion of 
the cusps from our considerations can be understood here. Although both and J| are 
well defined (zero) at the cusps, the left and right limits of , p and , z differ by 

V p' 2 +z' 2 \J p' 2 +z' 2 

an overall sign at the cusps (only —j == = is continuous at -0 = where it vanishes). Hence, 
e n and e^. are not well defined at the cusps. 

— * — * 

Because X is constrained to move in the surface, the differential dX, which is a vector 
with one-form coefficients, must lie in the tangent plane: 

— * — * — * — * — * 

dX = X^dip* +X >< f,d<f> = \\X^\\dip*e^ + ||X^||d0e^, 

= O^e^ + 9^ . (85) 
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The differential one-forms on the caustic surface are therefore 



, /l + C 2 (l-C 2 ) 



cos ^* 



sm- 



3^* 



dip* 



0<t> = pd> 

e n = o . 



a + - cos-^*(l + cos*0*) 



(86) 

(87) 
(88) 



The one-form 9^ vanishes only at the cusps. The line element, or the metric, ds 2 = dX-dX 
6i[>* ® @ip, + @<f> ® @<t> f° r the ring surface is 



ds = g^„d ip* + g^d <p= (p' + z' )d ip* + p d <\> 



-P 



i + C 2 , (i-C 2 ) , 

— 7, ^ 7> cos ^* 



sm 



d 2 ip* + 



P 



1 2 



a + - cos-^*(l + cosV'* 



d l 



(89) 



The area A of the caustic surface can be calculated by integrating the area element, pro- 
portional to the square root of the determinant of the above metric, over the whole surface. 
For the triaxial caustic rings (( — 1), we find 

i-2-k r-2n 



A = J dip J dip* det g~ij = 8irp + ^ 



(90) 



The length i of the tricusp cross-section can be calculated by integrating the square root of 
the line element with a constant 0, along the cross-section. For the triaxial caustic rings, 
we find 

r2n 

£= / dip*Vds 2 = 4p. (91) 
Jo 

Each vector field ei(ip*) is a differentiate map, and the differential dii at point P is a 
linear map. Thus, for each P, we can write 



(92) 



Since is a differentiable vector field, the lo^ are differential one-forms. They are called the 
connection forms in the moving frame ej. For the caustic ring, we have 



de 4 
de r . 



Vp' 



p ,2 + z ,2 

p'd(p 



p'd<p 

z'dcP _ _ _ 



VP' 2 + z 



(p'z" - z'p")dip* _ 
p' 2 + z' 2 6n 



Vp' 2 + z 



12 



(93) 
(94) 
(95) 
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Not all of the connection forms are independent. Differentiating the ortho normality condi- 
tion 81, and using Eq. 92, one finds that the connection forms are antisymmetric. Therefore, 

p'dcj) p'O^ (sin ip* + sm2ip*)d<j> 



(96) 
(98) 



Vp' 2 + z' 2 p^p' 2 


+ z' 2 v 


^[l+e+a-c^cos^*] 


sin^ 


{p'z" - z'p")d^ 


(p'z" - z'p" 






p' 2 + z' 2 


(p' 2 + z' 2 ) 


I ~" l + C 2 + (l-C 2 )cos^ 


z'd(f) z'e, 


b 


C(cos-?/>* — cos 2ip^d(f) 





^« ^ + Z n py/f/ 2 + z n ^ 2 [1 + C 2 + (1 - C 2 ) cos^*] sin ^ 

The connection one-forms uj^ and cj^„ are not well defined at the cusps — except for 
at ip — 0, where it vanishes. Although finite, their left and right limits differ by an overall 
sign at the cusps. cu^ n , on the other hand, is continuous everywhere. 

Equations 85 and 92, that define 9i and lo^ via dX and dii, describe how the moving frame 
varies as we move along a curve on the surface. The identities d(dX) = and d(dei) = 
(Poincare Lemma) imply 

= d(dX) = d(6 i e, i ) = dd^ - 0* A dii = (dO* - cu i:j A 9^ (99) 
= d(dii) = d(ujij(ij) = du^ej — lo^ A dij = (du>ij — oj^ A O0kj)&j ■ (100) 

The last equations in 99 and 100 that the forms and connection coefficients satisfy, 

dOi = ujij A 9j , dujij = Uik A uj kj , (101) 

are called the Cartan's first and second structure equations, respectively. Since there is only 
one linearly independent two-form on a two dimensional surface, we have 

dw^ = U^n A Unj, = —KgO^ A 6*0 , (102) 

where the scalar k g is the Gaussian curvature, and is independent of the choice of e^. 
Similarly, the two form 

#V A u (t>n ~ @4> A to^n = —2H9^p A Qfy , (103) 

defines the scalar H, called the mean curvature of the surface. For the caustic ring Eqs. 
102-103 yield 

z' (p'z" -z'p") 2C 2 [l + C 2 + (l-C 2 )cos^]- 2 

i^g\y*) p ( p i2 + z ,2y p(l + 2cos^)[a+f cos^(l + cos^)] ' 1 ' 
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and 



1 fp'z" -z'p" 



+ 



1 



Sign(cos^* — cos 2-?/^) 



V^[1 + C 2 +(1-C 2 )C0S^]3 



sin 



[f + cos^(i + cosvgyc 2 + cot 2 m 



(105) 



respectively. The Gaussian curvature at the outer cusp k g (0) = — However, k g is not 
well defined at the two non-planer cusps. It diverges with differing signs in the left and right 
limits at each of them. The mean curvature H is also divergent at the cusps. 

The connection one-forms ujp n and oo^ n are linear combinations of 6^ and 9^, for two 
dimensional surfaces, in general. Because of the relation d6 n = = u) ni p A 6^ + uj n p A 6^ we 
have symmetry in the expansion coefficients, hence 



lu, \ 

Using Eqs. 102, 103 and 106, we find 



\°<t>) 



(106) 



h^hcfKf, — h 



■>p(l> 



det (hi 



H = ^(V/> + h H ) = itrace (h i:j ) . 



(107) 
(108) 



Characteristic roots of the symmetric matrix (hy) are called the principal curvatures Ki of 
the surface. (Geometrically, they are the inverse "radii" along the principal directions of the 
surface, see the Appendix). We consequently have 

2 



2HKi + kq = , 



where 



Hence 



Kip K(f> , H = -(/ty, + K^) . 



K t = H±^H 2 - K G 

The matrix (/i^) is diagonal for the caustic rings, thus 



<2C, 



sin 



3-0* 



(p' 2 + z' 2 )-2 p[l + C 2 + (l-C 2 )cos^.]5 

z' ( Sign(cos-(/>* — cos 2^) 



PVP' 2 + z' 2 [a + | cos^,(l + cosV'*)] VC 2 + cot 



2 0* 



(109) 

(110) 
(111) 

(112) 

(113) 
(114) 
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The principal curvature is always positive. Its power series expansions 

Or 1 , C(2 - , ^ 



k*M = ^.(2tt - tf) = ^— + ^ v * ' + Oif) , (115) 



and 



* 3^' * a 3 ' 3p(l+3C 2 )4 p(l+3C 2 )l p(l+3C 2 )5 
show that /c^, diverges as 4 at the cusps, k^, on the other hand, is negative in the regions 
< ip* < if an d if < V'* < 27T, and positive in the region < ^ < 4f . As can be seen 
from the expansion 

= k,(2tt - 0) = -^^y + 0(^ 3 ) , (117) 

K(f, vanishes at the outer cusp. Although discontinuous, it remains finite at the two non- 
planer cusps: 

K ,(| T „) = K ^ ± „) = * 8 ^ + W + (#) . (H8) 

3 3 (8a-p)Vl + 3C 2 (8a -p)(l + 3C 2 ) 5 

Finite left and right limits at the cusps, differ by a sign. 

The (normal) curvature K n at a point on a surface in the direction of a tangent line L is 
calculated by using the Euler's Theorem 

/€„(?/>*, a;) = ^(■0*) cos 2 + K$(ip*) sin 2 a; , (H9) 

where is the angle between L and a tangent line in the principal direction associated with 
k^. As is shown in the Appendix, the inverses of and are respectively the curvature 
radii and along the principal directions on the surface. In the next section, we study 
the gravitational lensing effects of the dark matter caustic rings. Equations 114 and 119 are 
used to find the curvature radii along the directions of the line of sights. 



V. GRAVITATIONAL LENSING BY CAUSTIC RINGS 

The Caustic Ring Model of galactic halos precisely predict the density (Sect. Ill) and the 
geometry (Sect. IV) of the CDM distribution in the caustic neighborhood. Gravitational 
lensing may be a useful tool to test these predictions. Dark matter caustics have calculable 
lensing signatures [8-11]. In Ref. [9], we derived the lensing equations for the outer and 
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ring caustics, and estimated the image magnification at a number of sample locations of 
the line of sight on the caustic surface. The gravitational lensing effects of a caustic surface 
are largest when the line of sight is near tangent to the surface, because the contrast in 
column density, defined in Eq. 124, is largest there [9]. The effects, as can be seen in Eq. 
137, are proportional to (i) the fold coefficient A at the point where the line of sight is near 
tangent to the caustic surface, (ii) the square root of the curvature radius R of the surface 
along the direction associated with the line of sight, and inversely proportional to (iii) the 
critical surface density E c , defined in Eq. 125. Our focus is upon the caustic rings in this 
paper. Because, the fold coefficient increases proportional to (Eq. 66 and 70) as one 
approaches to the cusps of the rings (where ip* = \ ± ± , the effects are larger 
near the cusps. Because the curvature radius is largest near the outer cusp in particular, 
where R$ ~ -i? -1 (Eq. 132), the effects are largest there. Near the non-planer cusps, on the 
other hand, ~const (Eqs. 133), hence the effects are not as large there. In this section, 
to obtain the magnification for a line of sight tangential at an arbitrary point — except at 
the cusps where the lensing effects are infinite if 5v — 0, and our formulation fails — near 
the caustic surface, we use the equations for the fold coefficient and the curvature radius 
along the line of sight, obtained in Sects. Ill and IV respectively, in the lensing equations 
of Ref. [9]. We estimate the magnification by the rings at cosmological distances, where the 
critical surface density S c is minimized (hence the lensing effects are maximized), and also 
by the nearby fifth caustic ring of our own galaxy, at several sample locations of the line of 
sight in the caustic neighborhood, near the cusps. Like Ref. [9], we consider only the line 
of sights that are parallel to the galactic plane of the caustic ring and near tangent to the 
surface at a given ■?/>*. Thus, all the line of sights lie in the plane z = z(ip*) in this paper. 
Unlike Ref. [9], we pick ip* near the cusps, to estimate the lensing effects. 

We consider three cases of gravitational lensing by a caustic ring surface. In all the cases 
considered the line of sight is near tangent to the surface where a simple fold catastrophe 
is located. The three cases are distinguished by the curvature of the surface at the tangent 
point in the direction of the line of sight. In the first case, the line of sight is near tangent 
to the caustic surface that curves toward the side with two extra flows; see Fig. 3. We call 
such a surface "concave." 

In the second case, the surface is "convex," (i.e., it curves away from the side with two 
extra flows); see Fig. 4. 
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FIG. 3: Lensing by a concave fold. The arc is the intersection of the caustic surface with the plane 
containing the normal (z) to the surface and the line of sight (y). The shaded area indicates the 
side with the two extra flows. 



A 



3> 



FIG. 4: Lensing by a convex fold. Same as Fig. 3 except that now the caustic surface curves away 
from the side with two extra flows. 

In the third case, the caustic surface has zero curvature at the tangent point (the radius 
of curvature is infinite), but the tangent line is entirely outside the side with two extra flows. 

Although we exclude the cusp locations (where the lensing effects diverge if Sv = 0) from 
our formulation, we pick sample points that are near (but sufficiently far from) the cusps. In 
the presence of finite velocity dispersion, we constrain the expected effects at the smoothed 
cusps, considering the smearing out of the caustic surfaces of the axion and WIMP flows. 
To find an upper bound, as explained in Sect. Ill A, we first set the velocity dispersion of 
the CDM flow to zero, and then estimate the effects at the points on the surface that are 
respectively A*0* = and A*0* = ^ radian away from the cusps of the axion and WIMP 
caustic rings. Similarly, we estimate the lower limits at the points that are A^* = radian 
away from the locations that cups would occur if velocity dispersion were zero. 

Let us summarize the gravitational lensing phenomenon briefly. In linear approximation, 
the deflection angle 6 of a light ray due to a gravitational field is given by 

= v|yVdy, (120) 

where V(x,y,z) is the Newtonian potential. We choose the y-axis in the direction of prop- 
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agation of light. Geometrically, 9 is related to the angular shift £j — £ s on the sky of the 
apparent direction of a source: 

0(&) = §^(d -S), (121) 
where -Ds and -D^s are the distances of the source to the observer and to the lens respectively. 

— * — * 

£s is the angular position of the source in the absence of the lens while is the angular 
position of the image with the lens present. The angles carry components in the x and z 

— » — # 

directions: 6 = (0 X ,8 Z ), £ = (£ x , etc. Unless otherwise stated, we mean by a vector, 
a quantity with components in the x and z directions. We have x = (x,z) = -Dl£j- It is 
convenient to introduce a 2D potential so that 9 = -§f^^^(^,i) where = D L V, 

and D L is the distance of the observer to the lens. Then, Eq. 121 becomes 

£ = & + V^(5)- (122) 

It gives the map £s(£j) from the image plane to the source plane. The inverse map may 
be one to one, or one to many. In the latter case, there are multiple images and infinite 
magnification when a pair of images merge. The potential ip obeys the Poisson equation: 



Vl^^f^|. (1*.) 



where £(£z" x , £/ z ) is the column density (i.e., the integral of the volume density along the 
line of sight): 

= J dyd(D L ^ Ix ,y,D L ^ Iz ) , (124) 
and S c is the critical surface density 

s -«Sk^ o347g/cm2 (^ G p c )' < i25 > 

A uniform sheet of density S c focuses radiation from the source to the observer. In general, 
the shift in image position, A£ = £/ — £5, is the gradient of a potential whose 2D Laplacian 
is the column density. For an arbitrary mass distribution, the procedure for calculating the 
shift involves two steps. First, the matter density is integrated along the line of sight to 
obtain the column density. Second, the potential is obtained by convoluting the column 
density with the 2D Green's function. In our applications, however, the caustic has contrast 
in only one of the dimensions transverse to the line of sight. The procedure can be simplified 
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by expressing the shift directly as an integral over the parameter space of the dark matter 
flow forming the caustic. In this paper, we use the equations we derived in [9] to estimate 
the lensing effects in Sects. VA-C. 

Gravitational lensing produces the map £s(£j), given in Eq. 122, of an object surface 
onto an image surface. The inverse map may be one to one, or one to many. In the latter 
case, there are multiple images and infinite magnification when a pair of images merge. The 
image structure, distortion, and magnification are given by the Jacobian matrix of the map 
from image to source: 

K a = g^- = 5 a - i>a > ( 126 ) 

where vpij = ■ Because gravitational lensing does not change surface brightness, the 

magnification M is the ratio of image area to source area. Therefore, 

M = ■ (127) 



To first order, for ipij <C 1, 



E 



M = 1 + V 2 ^ = 1 + 2— , (128) 

where we used Eq. 123 in the last equality. In the convex (Sect. VB) and zero curvature 
(Sect. VC) cases considered, a point source can have multiple images [9]. In those cases, 
when two images merge, the Jacobian of the map vanishes and the magnification diverges. As 
we will see in Sects. VA-C, the magnification is larger near the cusps. We find that, for the 
line of sights considered, the magnification is proportional to near the non-planer cusps 
where the case is either concave or convex, on the other hand, it is proportional to near 
the outer cusp, where the case is concave (recall that $ — > at the cusps). Unfortunately, 
the gravitational lensing due to the fifth caustic ring of the Milky Way, which is only 55 pc 
away from us, turns out to be too weak to be observed — even near the cusps — with current 
instruments (Sects. VA-C). To obtain the largest lensing effects, we wish to minimize S c , 
given in Eq. 125. For fixed Ds, the minimum occurs when the lens is situated half-way 
between the source and the observer. Also, D s should be as large as possible. In our 
estimates, to get large effects, we assume that the source is at cosmological distances (e.g., 
2D L = 2D LS = D s = lGpc, in which case S c = 1.39 g/cm 2 ). We find in Sects. VA-B that, 
near the cusps of the caustic rings at cosmological distances, the effects for point-like sources 
are quite promising. The images of extended sources may also show distortions that can be 
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unambiguously attributed to the dark matter caustics [8, 9]. Observation of the calculated 
lensing signatures would give strong evidence for caustics and CDM in galactic halos. 

As is mentioned earlier, the lensing by a caustic ring surface also depends on the square 
root of the curvature radius R along the line of sight L. The principal curvature radii R^ 
and R<p are the inverses of the principal curvatures and k^, respectively. As is shown in 
the Appendix, is the radius along the cross-sectional plane of the caustic ring, and, R^ is 
the radius along the direction perpendicular to the cross-sectional plane. Therefore, R can 
be obtained from Eq. 119 as 

1 cos 2 cj sin 2 cj 

R = -R- + -R-> (129) 
where oo is the angle between the line of sight and the direction associated with R^. In this 
paper, we adopt the convention that R is negative (positive) if, along the line of sight, the 
surface curves toward (away from) the side with two extra flows. If R is negative, the surface 
is called "concave." If R is positive, the surface is called "convex." R^ > everywhere, 
except at the cusps where it vanishes. R<f>, on the other hand, changes sign at the non-planer 
cusps. In the regions where < ip* < y and -y < < 2n, we have R^ < 0. Between these 
regions, where y < < 4f , we have R<p > 0. The behavior of R near the cusps, where 
we expect large lensing effects to occur, can be obtained by expanding R^ and R^ as power 
series in the neighborhood of the cusps. Expansions around ip — 0, y, and 4^ show that R^ 
vanishes at all the cusps linearly: 

itV(tf) = R42n - tf) = M _ 9P(2 + 0(tf 5 ) (130) 



i?^, on the other hand, diverges as at the outer cusp: 

R M = - J) = J±±I^1 + [g - 3C)(a + P) ± + Q(#) , (132) 

whereas it remains finite (although discontinuous) at the other two cusps: 



T *) = ± 0) = T g^iM±g _ (»»-»>>« + O(^) . (133) 

Thus, for any ^ in the regions where < ip* < |7r, or |7r < ^* < 2vr, there is a pair of lines 
of sight for which the curvature vanishes. They are at angles: 



uj = ± arctan J — (134) 

V Rip 
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relative to the cross-sectional plane. Gravitational lensing by a fold of zero curvature is 
discussed in Section VC. 

A. Lensing by a Concave Fold of a Caustic Ring 

When the line of sight is near tangent to a caustic surface which curves towards the side 
with two extra flows, we define the fold surface as concave (Fig. 3). To estimate the lensing 
by concave folds, in Ref. [9], we considered only the outer caustics which are topological 
spheres surrounding the galaxies. The outer caustics are concave everywhere on the surface. 
Here, in this paper, we study the gravitational lensing by the concave folds of the caustic 
rings. If the line of sight is chosen along the principal direction associated with R^, two 
thirds of a caustic ring surface, where < ■?/>* < an d < ip* < 2n provides a cross- 
section where the fold is concave; see Fig. 5. The convex and zero curvature cases are 
discussed in Sects. VB-C. 



A 



B 




v 



FIG. 5: Lensing by a concave fold of the caustic ring for a line of sight near tangent to the surface 
at a point ij)* close to the outer cusp, where (p, z) = (po, 0). The line of sight lies in the z = z(ip*)- 
plane. We define x = p — po- A) Side view in the direction of the line of sight. The latter is 
represented by the dot near x = z = 0. B) Top view. The solid curve is the location of the concave 
fold in the z = z(^*)-plane, and the dotted curve is the location of the outer cusps of the ring. 



38 



Image shift and magnification for a concave fold are given [9] as 



A£ = 6-e5 = 776©(-6) (135) 
M = § i = l + ^6(-6) + 0(r ? 2 ), (136) 



where 



2nAJ2 R 



V = * • (137) 



When the line of sight of a moving source crosses the surface of a simple concave fold, 
the component of its apparent velocity perpendicular to the fold changes abruptly. Also, a 
discontinuity occurs in the magnification of the image. Both effects are of order rj. Here, 
using Eqs. 66, 70, 132 and 133, we see that, for the concave case, the maximum effects 
should be expected near the cusps. Near the outer cusp rj ~ where ip* = $ or 2ir — d. 
Near the non-planer cusps rj ~ t9~^ where ip* = ^ — *& and ip* = 4^ + *&. Inserting the fold 
coefficient A for the rings, given in Eq. 61, and R = R^ip*) = /tl 1 ^*) given in Eq. 114, 
into Eq. 137, we obtain 

= 4^ 2 M cosa(^) 
' Y, c dttdT b^/Ep W *> ' 1 ; 



(139) 



where 

| esc ip* I 



^ [1 + ^ cosV*(l + cos^*)]|(l + 2cosV*) tan^j 
Notice that, ( dependencies of A, and cancel each other in Eq. 137, hence rj in Eq. 138 
is independent of (. Combining Eqs. 59 and 138 gives rj n: for the line of sights that are 
parallel to the galactic plane of the caustic ring under consideration and near tangent to its 
surface at any as 

*<«.) = ^U v f^^SM . (140) 

It was shown in reference [6] that b n and v n are of the same order of magnitude. Moreover, 
Sikivie [13] interpreted the ten rises in the rotation curve of the Milky Way as the effect of 
caustic rings. In that case, the widths p n of caustic rings are determined from the observed 
widths of the rises. Typically, one finds p n ~ 0.1 a n . Using this, v n ~ b n , and cos a(ip*) — 1) 
Eq. 140 yields 

{„„ : „ = 1,2, ...} ~ (7,6,6,5,5, ...) . 10"U(«^ g)(A)(_^_) . (141) 

Comparing the estimates Eq. 141 with the estimates of rj n for the outer caustics considered 
in [9], where the simple fold is always concave, we see that the lensing effects of the caustic 
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rings are about 10 G n {jp*) times larger than the lensing effects of the outer caustics. We may 
expand the function Q n (^) in Eq. 140, around ip* = 0, ^, and 4^. Near the outer cusp 

J2-&- 1 (5 + 14^)^ [23+4^(9 + 37^)]tf 3 

V ; V ; V3(l + ^)i 12^(1 + ^)1 320^(1 + ^)5 v ; v ; 

For i? <C 1, the first term dominates, hence increases proportional to i? -1 , as i? 

decreases. For p n = 0.1 a n , Q{"&) ^ 0.778 Near the non-planer cusps (where the cases 
are concave) we have 

2tt 4tt 4iH tf5 (3-^2n)^f 

a „ ( _ _ = a „(_ + „, = + + + 0( „ 5) . (143 , 

For small the first term dominates. Taking p n = 0.1 a n , we find Q{^f — i?) — 0.624 • 
The expansions for Q n {^f + #) and G n {^f — will be used in Sect. VB, where the cases 
are convex. We investigate the gravitational lensing properties near the cusps of axion and 
WIMP caustic rings in the next two sections. As in Sect. Ill A, to estimate the upper 
bounds at the smoothed cusps (when 5v ^ 0), we pick d = and $ = ^ (implied by the 
minimum primordial 5v of the particles) for the axion and WIMP caustic rings, respectively. 
The lower bounds are estimated choosing ■& = (implied by the triangular feature [13] in 
the IRAS map of the Milky Way galactic plane). 



1. Lensing by a concave fold near the cusps of axion caustic rings 

At the outer cusp of an axion caustic ring, taking p n — 0.1 a n , we constrain the range of 
function Q n between (/(i^) ~ 1.9 and G(±-^!—) ~ 1858.5 (for the fifth ring of the Milky 
Way, we have G 5 (±y^) ~ 2, and Q 5 (±j^) ~ 1933.7). Assuming that the line of sights 
are parallel to the galactic plane of the caustic ring, and near tangent to the surface at the 
points ■?/>* = ±7^)0 an d i ) * — ±75 respectively; see Fig. 5, we find r](±j^) ~ 2.91 ■ 10~ 2 and 
r](±j^) ~ 27.88 for the sources at cosmological distances (e.g. 2D L = 2D LS = D s = 1 
Gpc). Thus, depending on the magnitude of the effective velocity dispersion of the axion 
flow in a galactic halo, 77 may range between 0.03 and 28 at the outer cusp of an axion caustic 
ring. This means that the magnification of a point-like source can be between 3% and 2800%. 
The upper bound, however, is estimated at a point whose distance to the location that the 
cusp would occur if 8v were zero, is about the minimum smearing out distance 5x a due to 
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the primordial velocity dispersion of the flow (see Sect. Ill A). Because the effective velocity- 
dispersion of the flow is expected to be larger than its primordial velocity dispersion, caustics 
are likely to be smeared out more than the size 5x a estimated in Eq. 5. To make another 
estimation at a point further than ^ = from the outer cusp, let us choose a sample point 
at ip* = j^q, where Q{-^ ) — 24.79. For the sources at cosmological distances, we obtain 
77(1^0) — 0.37, which implies about 37% magnification. This is about 250 times larger than 
the effect estimated (of order 0.15%) for the concave folds [9] considering the outer caustics 
at cosmological distances. The 37% magnification is about 37 times larger than the largest 
effect we estimated (of order 1%) in Ref. [9], considering a line of sight near tangent to 
a cosmological caustic ring at (p,z) = (a, 0), where the fold is convex. At this location, 
however, the lensing effect of the convex fold is minimum, because, the fold coefficient A 
and the curvature radius R<p are minimum there. Both the A and increase monotonically, 
as the line of sight approaches to the cusps. At the non-planer cusps of a cosmological axion 
caustic ring, we have r/(f -^) = rj (f + ^) ~ 0.016, and vif-^oo) =*Kf +tB5o) - 0M - 
This implies that the magnification at the non-planer cusps of the axion caustic rings may 
range between %2 and %46. At sample points that are ^0 radian away from the non-planer 
cusps (i.e. choosing d = ^ in Eq. 143), we find G = G(*$j&) ^ 3.54. Therefore 
^(^f) = ^(Hlf) — 0-053. Hence, the magnification is about 5%, at the sample points near 
the non-planer cusps of the cosmological axion caustic rings, when the case is concave. 

Let us estimate the 77, for the fifth caustic ring of the Milky Way, which is the one closest 
to us and believed to be the most constrained by observation. We have a 5 = 8.31 kpc, 
b 5 = 657 km/s (if r > 0) or 516 km/s (if r < 0), p 5 = 0.134 kpc, q 5 = 0.2 kpc, f 5 = 0.02 
and t>5 = 480 km/s. The closest distance to the fifth ring is 55 pc. The tangential distance 
to the ring, Dl, is therefore 965 pc. Thus 

where is about 1.37 if r > 0, or about 1.27 if r < 0. Considering the cosmological 

sources for which D^s/Ds — 1, we find 775 (V**) ~ 1-7 ■ 10~ 7 (?5(-?/>*). Near the outer cusp 
^(t?) ~ 0.81 • -i?" 1 . Unfortunately, even at = (assuming the ring is axionic and the 
velocity dispersion of the flow is about the primordial value), we find that 775 ~ 3.33 • 10~ 4 
if t > 0, or r/ 5 ~ 3.10 • 10~ 4 if r < 0. Thus, the magnification and image distortion are 
negligible (of order %0.03). However, as noted in Refs. [8, 9], the images of extended sources 
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may be modified in recognizable ways. In particular, straight jets would be seen with an 
abrupt bend where their line of sight crosses a fold. If a jet makes an angle a with the 



At the outer cusp of a WIMP caustic ring, taking p n = 0.1 a n , we constrain the range 
of function Q n between Q(±y^) ~ 1.9 and Q(±j^) ~ 18.6 (for the fifth ring of the Milky 
Way, we have Q 5 (±-^) ~ 2 and C/ 5 (±^) ~ 19.3). Assuming that the line of sights are 
parallel to the galactic plane of the caustic ring, and tangent to the surface at the points 
V>* = ±YE and ^* = ± # respectively, we find ^(±7%) ^ 0.029 and ri(±£) ~ 0.28 for 
the sources at cosmological distances. Thus, depending on the magnitude of the effective 
velocity dispersion of the axion flow in a galactic halo, rj can be in the range between 0.03 
and 0.28 at the outer cusp of a WIMP caustic ring. This implies between 3% and 28% 
magnification near the outer cusp. At the non-planer cusps of a cosmological WIMP caustic 



ring, we have r](f - jj) = v(f + Ts) - °- 016 and v(f ~ #) = v(f + - °- 046 - This 



implies that the magnification at the non-planer cusps of WIMP caustic rings may range 
between %2 and %5 when the fold is concave. If the nearby ring is a WIMP caustic, its 
lensing effects near the cusps are negligible (at most, we find ^(i^) ~ 3.3 • 10~ 6 , and 
775 (^f — |g) = 775 + j^) ~ 5.2 • 10~ 7 ). In the next section, we consider the gravitational 
lensing effects of caustics rings when the line of sight is tangent to the surface where a convex 
fold catastrophe is located. 



By definition, a convex fold is curved in the direction opposite to the side with two extra 
flows (Fig. 4) and radius of curvature of the surface along the line of sight is positive. For 
a caustic ring, in the region where ^ < < 4^, > 0, hence the simple folds associated 
with in this region are convex; see Fig. 4. We assume that the line of sight is parallel 
to the galactic plane of the caustic ring, and passes by the surface near tangent as in Fig. 4 
(hence the line of sight lies in the z = z(^)-p\ane). The ^/-dependent shift is given [9] as 



normal, it appears bent [9] by an angle 5 = \r\ sin 2a. 



2. Lensing by a concave fold near the cusps of WIMP caustic rings 



B. Lensing by a Convex Fold of a Caustic Ring 




(145) 
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where the length scale L is a cutoff beyond which our description of the flow is invalid. 
The function rj is given in Eq. 140. The ^/-independent part of the shift, — j£^ R , is 
uninteresting, and therefore subtracted out in Eq. 145. The magnification is given [9] as 





-1 


dii 





1 + v_ ln [R,D L \o 



(146) 



7T V 2L 2 

The cutoff L has an effect on the magnification and elongation of the image in the direction 
normal to the caustic surface, but that effect is independent. L has a global effect on 
the image, as opposed to an effect localized near = 0. In particular, when the source is 
exactly behind the caustic (£g = 0), the images are at £/ = — £ c , 0, and +£ c , with 

£ C =Tr-^exp(-- + l). (147) 

A point-like source has a single image sufficiently far from the caustic, say at £n > 0. When 
the line of sight approaches the caustic surface tangent point, two new images appear on top 
of each other at £/ 2 = £13 = — £c/e- At that moment, the magnification at £/ 2 is infinite, and 
£,s = V^c/en. As the source crosses the caustic, £/ 2 moves toward £n and finally merges with 
it. When £n = £/ 2 — +£c/ e > the magnification diverges again. After that, only the image at 
£/3 remains. The image of an extended object is modified [9] in the direction perpendicular 
to the convex fold by the relative amount 



where 



X_i = _?Zm(M , (148) 

71 \ZdJ 



21 2 

(149) 

The image is stretched for £/ < £ d , and compressed for > £ d . The line of sights we 
consider in this section are near tangent to the caustic surface at sample points ^ = ^ + -d 
(ip* = 4f — $)j where 1? are small, and lie in the z = z{^f + i?)-plane (z = z(^f — $)-plane). 
For such a line of sight, L 2 is of order \^3(8a — p)p$ 3 /IG. Therefore, using Eqs. 133 and 
149, we obtain £ d ~ 3(p$ 3 /^l + K 2 D L . Taking D L = 0.5 Gpc, p = 100 pc and ( = 1, we 
estimate C, d ~ 3- 10~ 7, $ 3 at the sample locations of cosmological caustics. For the fifth caustic 
ring of the Milky Way, we find ^5 ~ 0.2 • $ 3 . Because p/Di is the transverse angular size of 
the caustic ring, our description certainly fails for £/ > p/D L . We have p/D L ~ 2 • 10~ 7 for 
the rings at cosmological distances, whereas p$/D L5 ~ 0.14 for the fifth ring of the Milky 
Way. 
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The estimates for 77 (Eq. 137) in Eqs. 145 and 146 are the same as in Eq. 141, except 
that, now, G(ip#) has to be evaluated in the region where ^ < ip* < Note that at ip* — n, 
i.e at (p, z) = (a, 0), for p n — 0.1 a n we have Q n (7r) = -75, and hence, rj n {n) are equal to the 
estimates obtained (of order 10 -2 ) for the convex case in Ref. [9] (denoted by r]' n ). However, 
near the non-planer cusps, we have 

&( _ + 1) = w _ _ v = _ + + o,,,) . (150 ) 

Because Q n diverge as $~2 near the non-planer cusps, r] n (ip^) are larger than i] n (7i) near the 
cusps. Depending on the magnitude of the effective velocity dispersion of the CDM flow in a 
galactic halo, the range of r\ can be constrained at the smoothed non-planer cusps. For the 
axion caustic rings, taking p n = 0.1 a n , we find 77 + yg^j) = 0.46 and rj(^- + yy) = 0.03, 
for the rings at cosmological distances. (Due to the reflection symmetry of the caustic rings, 
the same results hold at tp* = — yg^j and ty* = — j% respectively). For the WIMP 
caustic rings, taking p n — 0.1 a n , we constrain the range of 77 at the smoothed non-planer 
cusps between 77 + yr) ~ 0.05 and + yy) ~ 0.03, for the sources at cosmological 
distances. (Due to the reflection symmetry of the caustic rings, the same results hold for 
ip* = 4y — ^ and ip* = ^ — yy). Unfortunately, although rj increases as the point at which 
the line of sight is near tangent to the convex fold approaches to the non-planer cusps, the 
magnification (Eq. 148) is not. This is because, the line of sight stays close over shorter 
depths, as it approaches to the non-planer cusps. For an angular distance £/ ~ 10~ 9 , we 
find about 1% magnification at the sample point tp* = ^f - + for the cosmological caustic 
rings. Even at angular distances as small as ~ 10~ 10 , the magnification is about 3% at 
this point. The lensing effects by the convex fold of the nearby fifth ring are negligible (of 
order 10~ 6 ). Finally, in the next section, we consider the gravitational lensing by a fold with 
zero curvature. 



C. Lensing by a Zero Curvature Fold of a Caustic Ring 

We found in Eq. 134 that, in the regions < ip* < \ an d -f- < ^* < 2n a caustic ring 
has a pair of tangent lines along which the curvature vanishes. For such a line of sight, the 
image shift is given [9] as 

A£ = -e(-0)(-£o0 3 ) 1/4 , (151) 
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where 

£o = (9-89 A)4iL . (152) 

Here, U is a positive constant and has dimensions of (length) 3 [9]. The fold coefficient A is 
estimated in Eq. 63. The magnification is given [9] as: 

-i 

(153) 

Triple images occur when < £o- Unfortunately, for the zero curvature tangents of 
caustic rings, £o is very small. To have some order of estimates for let us consider the 
line of sights along the zero curvature tangents at tp* — ^ and ip* = of a cosmological 
ring surface (2D L = 2D LS = D s = Gpc). From Eq. 63, we find A(^) ~ 1.5 • 10~ 3 , gr L 

cm 2 kpc2 

and A(f=) ~ 5.1 • 1(T 4 gr 1 . Taking U = (kpc) 3 , one finds ^ 2.9 • 10" 14 and 

cm 2 kpc2 

£0(715) — 3.5 ■ 10~ 16 . Hence, the triple images cannot be resolved. At an angular distance as 
small as ~ 10~ 9 (10 -10 ) the magnification and image distortion are of order 2% (3%) for 
the line of sight near tangent at ■?/>* = 75 • For the line of sight near tangent at ^ = the 
effects are of order 6% for ~ 10~ 9 and 11% for ~ 10~ 10 . For the fifth caustic ring of 
the Milky Way the effects are negligible (of order 10" 5 for ~ 1CT 9 and ip* — ^). 

In the next section, we study the structural stability of the caustic rings, using the tools 
of the Catastrophe Theory. 



m=i- feK,) (-§ 
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VI. STRUCTURAL STABILITY OF CAUSTIC RINGS 

In this section, we analyze the dark matter caustic rings in the Catastrophe Theory [19] 
point of view. Catastrophe Theory studies phenomena distinguished by sudden changes in 
behavior, arising from smooth alterations in circumstances. It analyzes how the qualitative 
nature of equation solutions that describe the states of the phenomena, depends on the pa- 
rameters that appear in the equations. In general, the problem of determining the solutions 
(states), let alone analyzing how these solutions change, as the parameters that control the 
circumstances change, is a formidable task. Considerable amount of simplifications occur, 
however, if the equations can be derived as the gradient (with respect to the state variables) 
of a catastrophe function (or, potential function) $ of both the state variables and control 
parameters. In particular, Catastrophe Theory studies how the equilibria (or the critical 
points), V$ = 0, of the gradient systems change, as the control parameters change. 
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The catastrophe function of the triaxial caustic rings (( = 1), whose equilibrium points 
are given by the flow equations (Eqs. 12-13) near the caustic, can be given as 



^(XuX2,P,P~a,z) = - Xi 3 - X1X2 2 + V / PX2 2 - (p - a)xi ~ zXi ■ 



(154) 



In the Catastrophe Theory terminology, the x% are called the state variables, p, the "lon- 
gitude" of the caustic cross-section, and (p — a,z), the physical space coordinates on the 
transverse cross-sectional plane (Fig. 1), are called the control parameters of the $. The 
term | xi 3 ~ X1X2 2 , which is purely composed of the state variables, is called the germ, and 
the terms proportional to the control parameters are called the perturbations. We have seen 
that neither the existence of the caustic rings (Eq. 21), nor their lensing effects (Eq. 138) 
depend on (. (Both the function 77 and the vanishing of determinant D 2 are independent of 
(). The parameter ( determines the ratio of the latitudinal dimension q, and the longitudi- 
nal dimension p. Recall, however, that the flow equations near the caustic rings are derived 
[6] assuming that the caustics are tight, i.e. p <C a and q <C a. The fifth caustic ring of 
the Milky Way appears tight in the IRAS map [13], and ( 5 = = 115 in that case. 

When ( = 1 (£ 7^ 1), the three dual cusps of the caustic cross-section make up an equilateral 
(arbitrary) triangle. 

The potential given in Eq. 154 is symmetric under the discrete transformations (xi, X2) —> 
(Xi, — Xi) and (p,p — a,z) — > (p, p — a, —z). The equilibrium, V$ = 0, of the caustic ring 
potential indeed implies the flow equations (Eq. 12-13 with ( = 1) 

.9$ 



dxi 

•9X2 








1 2 2 
a + Xi - X 2 



2(v / P- Xi)Xi ■ 



(155) 
(156) 



If the system described by $ is in equilibrium (stable or unstable) at a particular point 
in the state space, the stability properties of the equilibrium can be determined from the 
Stability (Hessian) Matrix 



9$ 
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VP- Xi ) 



(157) 



This is, in fact, nothing but the Jacobian matrix we have obtained in Eq. 20. The equi- 
librium points, or the critical points, at which V$ = are called nonisolated, degenerate, 
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or non-Morse critical points. The critical points at which det $jj 7^ are called isolated, 
nondegenerate, or Morse critical points. The fourfold degenerate critical point of $ occurs 
when the stability matrix $^ vanishes identically. This implies xi — X2 = P = 0. Because 
all the critical points a priori satisfy Eqs. 155-156, we also have p = a, and z = at the 
fourfold degenerate critical point. Thus, the point at (p,p,z) = (0, a, 0) G R 3 in control 
parameter space, parameterizes the function 

$(Xi,X 2 ,0,0,0) = ^Xi-XiX2, (158) 

which has fourfold degenerate critical point at (xi: X2) — (0, 0) G R 2 , in state variable space. 

The twofold and threefold degenerate critical points are found by requiring one of the 
eigenvalues of the stability matrix to vanish. When this occurs, the determinant of $y is 
zero: 



£>2(xi,X2) = -4 



Xi --3- I +X2 - i 



. (159) 



So, the set of degenerate (Non-Morse) critical points of $ is the circle 

This means that, whenever a critical point lies on this circle, it is doubly or triply degenerate 
(fourfold if p — also). As we have seen in Sect. II, points of this circle in the state variable 
space, produce the caustic cross-section with longitude p in physical (p — a, z)-plane. Now, 
to study the stability properties of the caustics (degenerate critical points), we can use the 
two state variables xi an d X2, to give a parametric representation for the three control 
parameters in R 3 . The parametric representation for p(xi>X2) that can be obtained from 
Eq. 160 is 

P={^ A ) 2 ■ (161) 

The parametric representation for p and z are given in Eqs. 155-156. Instead of working 
with these parametric equations directly, to study the stability properties easily, we may 
take advantage of the scaling relations among the state and control parameters. As can be 
seen through Eqs. 154-157, if the state variables xi and X2 both scale like A, then all the 
control parameters p, p — a, and z scale like A 2 : 

(Xi, X2) -> A(xi, X2) =>(p,p- a, z) -> X 2 (p, p-a,z) . (162) 
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Moreover, because the longitudinal size of the caustic rings p > 0, it is sufficient to consider 
the two dimensional (p — a, z) cross-sections of the three dimensional control parameter 
space, for p = 0, and for any p > (recall that we study the tight caustics for which 
< p <C a, in this paper). Any other case can be considered via the scaling relations Eq. 

162. For p > 0, we determine how the circular set (Eq. 160) of degenerate critical points 
parameterizes the (p — a, z) cross-section by substituting a parameter representation (Eq. 
22) of this circle: 

Xi = ^(l±cosV>) » X2 = ^sin^, 

where ip e [0, 2n] is an angular variable, into Eqs. 155-156. We find 

p p 
p(ip) — a + - cos'0(cos'0 ± 1) , z(ip) = - sin-^(l =F costp) ■ (163) 

These are the caustic equations (Eqs. 23-24) with £ = 1. Thus, by substituting the paramet- 
ric equations (Eq. 22) for Xi an d X2 satisfying the two and three fold degeneracy condition 
(det <3>y = 0), into the equations for critical points (flow equations), we have found, in Eq. 

163, the doubly and triply degenerate critical points (folds and cusps respectively) of the 
catastrophe function <3>. The symmetry (xi, X2) (xi, — X2) an d (p, p—a, z) — * (p, p—a, —z), 
which was noted earlier, is realized as ip — > —ip in Eqs. 160 and 163. When p = 0, Eqs. 155- 
156 and 160 imply that the caustic cross-section reduces to the point (0, 0) in the (p — a, z)- 
plane. In this case, the caustic tube becomes an isolated caustic circle with (p,z) = (a,0). 
(It becomes a topological circle if the axial symmetry is not assumed). The plot of Eqs. 163, 
for a given p > 0, in (p — a, z)-plane is called the -D_4 catastrophe which we call tricusp; 
see Fig. 1. For any p > 0, the tricusp on the (p — a, z) -plane of the control parameter 
space, divides the plane into three open regions which describe catastrophe functions of 
three qualitatively different types. The qualitative properties of functions parameterized by 
points within any one region, however, are the same. These properties change, as we pass 
through the tricusp on the plane. Therefore, to determine the qualitative properties of the 
catastrophe functions, it is sufficient to consider convenient points in each of these regions. 
We restrict ourselves to the p — a axis. On this axis, one passes through the tricusp at 
p — a = and p = po = a + p. Hence, the three open regions are: p < a, a < p < po, and 
p > po; see Fig. 1. The equations 155-156, which determine the critical points of $, on this 
axis are: 

p-a-xl + xl = (164) 
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2( v / P-Xi)X2 = 0. (165) 

We infer, from Eq. 165 that the critical points must have either (i) X2 = 0, or (ii) xi — y/P- 
If (i) is the case, Eq. 164 yields Xi — ±VP — a ( out an d m flows) for p > a. If (ii) 
is the case, Eq. 164 yields X2 = ±\/Po ~ P (down and up flows) for p < p . Thus, for 
any given p > 0, along the p — a axis of the control parameter space of the catastrophe 
function $, there are two sets of critical points for p < a: (xi, X2) = (\/P, ±\/Po _ p)(down 
and up flows), four sets of critical points for a < p < p : (xi,X2) — (y/P, i^Po — P) and 
(.X11X2) = (i^Po — P, 0) (down, up, in and out flows); see Fig. 6, and two sets of critical 




P<Po a<P<P 



FIG. 6: The critical points of the catastrophe function for the caustic rings are shown in the (xi- 
X2)-plane as a function of the control parameter p for given ring radius a and caustic longitude p. 
In the left panel p < po, whereas on the right panel a < p < pq. The circular sets of degenerate 
critical points are located where the determinant of the stability matrix vanishes. The S and M 
denote the saddle and local minimum respectively. The signs of the eigenvalues along the principal 
directions at the isolated critical points are indicated by the thick arrows. 

points for p > p : (xi,X2) — (±VPo — P, 0) (in and out flows); see Fig. 7. 

Once the locations of the critical points are known, their stability properties can be 
determined from the Stability Matrix Eq. 157. In case (i), where (xi,X2) — (±-\/p — a > 0), 
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FIG. 7: The critical points of the catastrophe function for the caustic rings. Same as Fig. 6 except 
that now the control parameter p > pq. 



(for p > a), we have 



*ij = 2 



±y/p^d o 

y/P =F VP - a 



(166) 



Thus, in this case the xi an d the X2 axes are the principal directions with the eigenvalues 
A X1 = ±2 y/p — a and A X2 = 2( v /p=p y/p — a) respectively. The xi eigenvalue of the upper 
(lower) critical point A^ (A Xl ) is always positive (negative), except at p = a, where it 
vanishes. The xi eigenvalue of the upper critical point, A^ 2 = 2(^/p — y/ p — a) changes sign 
when p goes through p = a + p: A^ 2 is positive for a < p < p , whereas it is negative 
for p > pq. Thus, the upper critical point (xi,X2) = (y/p — a, 0) with the eigenvalues 
(A^,A^ 2 ) = 2(y/p — a, y/p — y/p — a) is a local minimum for a < p < po, whereas it is 
a saddle for p > p (A^ 2 vanishes when p = p ); see Figs. 6-7. The lower critical point 



(Xl,X2) 



-y/p — a, 0), on the other hand, has the Xi eigenvalue A^ = —y/p — a < 



(vanishes when p = a), and X2 eigenvalue A^ 2 = y/p + y/p — a > 0. Thus, the lower critical 
point is always a saddle for any p > a. Note that, the points p = a and p = p are the 
locations at which the critical points (xi,X2) = (±VP — a i 0) cross the circular degenerate 
critical set xi + x\ ~ y/PXi — (Eq. 160) on the state variable plane. We, therefore, expect 
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these locations to be special. 

The stability properties and principal axes of the critical points of case (ii), (xi,X2) — 
(y/p, ±VPo — p), where p < p , are determined similarly. The stability matrix Eq. 157, for 
this case is 



(167) 



y/P TVPO ~ P 

tVpo - P o 

The principal directions are no longer along the Xi an d X2 axes. So, they must be determined 
by matrix diagonalization. We find that, the eigenvalues Ai and A2 associated with the 
principal directions Vi and v 2 are the same for the upper and lower critical points: 

(K,K) = (OU = (y/p + y/p + 4(p -p),Jp-y/p + 4(p - p)) . (168) 

The eigendirections at the upper and lower critical points are 

/ 1 \ 



2y/p -p 



y/P+yJp+HPQ-p) 

2%/po-P 



(169) 



and 



\ 



V 



-\/P+y/p+Hpo-p) 



4 



2^/po-p 



\ 



V 



fP- y/p+ 4 (p»-p) 
2Vpo-p 



(170) 



A 



VI 



y/p + \/p + 4(p — p) > 0, whereas A 



V2 



A! 



respectively. The eigenvalues A^ 

y/p — Jp + 4(p — p) < (vanishes at p = p ). Thus, both of the critical points (xi, X2 



V2 



(y/p, ±VPo — P~), where p < p , are saddle points; see Fig. 6. Where p = po, the critical 
points cross the circular degenerate critical set in the state variable plane, and, A" 2 and \ l V2 
vanish. We expect this point to be special too. Before investigating the special locations 
at p = a, and at p = p , let us summarize what we have learned. In each of the regions 
divided by the triangular shape of the (p — a, z) cross-section of the control parameter space, 
different types of catastrophe functions $ are parameterized. The interior of the triangular 
region (the region where a < p < po 011 the p-axis) parameterizes the functions which 
have three saddles, and one local minimum. The saddles exist outside the circular set of 
degenerate (Non-Morse) critical points, and the local minimum exists inside the circular set 
(Figs. 6-7). The region outside the triangular domain (where p < a and p > p on the 
p-axis) parameterizes the functions with two saddles, both of which lie outside the circular 
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set of degenerate critical points of the state variable plane. Let us, now, investigate what 
happens at the special locations p = a and p = p , as p increases along the line z = for a 
given p > 0. For p < a, there is a pair of saddle points at (xi, X2) = (y/p, ±\/Po — p), each 



of which has the eigenvalues X VuV2 = A /p ± \Jp + 4(p — p) in the corresponding principal 
directions. As p increases through a, a fold catastrophe occurs at (X11X2) = (0,0) (when 
p = a). This is because, once p > a, a saddle comes into existence outside the circular set, at 
(Xi> X2) — (— \J P — a i 0) with (A Xl , A X2 ) = 1{—\Jp — a, y/p + \fp — a), and a local minimum 
comes into existence inside the circular set, at (xi,X2) — (y/p ~ a i 0) with (A Xl ,A X2 ) = 
2(y/p — a, y/p — y/p — a). The new critical points at (xi, X2) = (±VP — a i 0) approach each 
other as p approaches a from above and become degenerate at (xi, X2) — (0, 0) when p = a, 
and disappear for p < a. Meanwhile, the eigenvalues in the X2 direction approach to ^/p, 
and the eigenvalues in the Xi direction decrease monotonically, and vanish at p = a. This 
can also be seen from the fact that, at the doubly degenerate critical point (xi, X2) = (0, 0) 
(where p = a on the p-axis), the catastrophe function becomes 

$ = \xl - xixl + VpxI , (171) 

and the stability matrix takes the form 

(° \ 

®a = ■ 172 

Thus, the eigenvalue in the xi direction is zero, whereas the eigenvalue in the X2 direction 
is nonzero and positive. There exists a catastrophe at p = a, because, the number of 
isolated critical points in the (xi, X2)-plane changes as p passes through a, where det $jj = 
0. Moreover, the point p = a on the p-axis is expected to be a location for a fold (A 2 ) 
catastrophe, because two critical points are involved there. The rigorous way to prove this, 
is to check if the form of the catastrophe function Eq. 171 can be reduced to the canonical 
form for an A 2 catastrophe, by making a smooth change of variables at p = a. Since we 
expect the catastrophe germ is of type A 2) which is C2X2 + c iXi, an d from Eq. 154 we have 
yjp as the coefficient of xh ^ e canonical form at p = a that we seek is y/px'i + c iXi- We 
introduce the following non-linear transformation 

Xi = Xi + Ai = xi + A 2oxl + A nXiX2 + ^02X2 

X2 = X2 + A 2 = X2 + B 20 xl + £11X1X2 + £02X2 ' ( 173 ' ) 
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and see if the coefficients A^, and can be chosen to satisfy 

VpxI - Xixl + 7^x1 = x/px'2 + c iXi 3 

= VPixl + 2% 2 A 2 + Aj) + Cl (xl + 3 X ?A! + 3 X iA? + Af) . (f 74) 

The equality is already satisfied for the terms of degree two. For the terms of degree three, 
we have 

-Xixl + \xl = 2v^X 2 A 2 + c lX l = 2^(^20X1X2 + #11X1X2 + #02X2) + CiX? • (175) 
The solution is 

#20 = 0, 5 " = 2^> ^02 = 0, ci = -. (176) 

On the left hand side of Eq. 174, the coefficients for the terms of degree four and higher are 
all zero. The corresponding terms on the right hand side of Eq. 174, however, depend on 
the disposable coefficients Aij and Bij (for i + j > 2). Because, for a given degree, there are 
more disposable coefficients on the right hand side of Eq. 174, than the zero coefficients of 
the left hand side, it is possible to choose and B^ so that all the monomials of degree 
greater than or equal to four, have vanishing coefficients. Thus, at p = a, the catastrophe 
function can be put into the form 

$ = VPX'l + \x'? • (177) 

By rescaling the control parameters p^x'2 = V an d ^~^Xi = x i we obtain the canonical form 
of A 2 catastrophe: $ = y 2 + x 3 . This proves that p = a is a fold location. 

We now turn our attention to the other interesting location on the p-axis, that is, to 
the point p = p . As p increases further toward po, the new saddle on the xi ax is at 
(Xi,X2) = {-y/p- a,0) with (A X1 ,A X2 ) = 2(- v /p-a, v /p - yj p = a ) moves away to the 
west. The local minimum (xi, X2) = (y/p — a, 0) with (A X1 , A X2 ) = 2(^p — a, y/p — y/p — a) 
moves across the degenerate critical circle, toward the point (xi,X2) = (y/P,ty- The two 
original saddles at (xi,X2) = (y/p, ±\/Po ~ p), with A^ = y/p ± yp + 4(p - p), move 
toward the same point (y/p, 0) along the line xi — yfP- As the critical points approach to 
(y/p, 0), the two principal directions of the two saddles rotate toward the xi an d X2 axes, and 
the eigenvalues all approach to (2 y/p, 0). As p passes through p , the two original saddles 
and the local minimum collide at (xi,X2) = (y/p,Q), and form a triply degenerate critical 
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point. By expanding the catastrophe function Eq. 154 for z = and p = p about the three 
fold degenerate critical point (xi, X2) = (\/Vi 0) using xi = (xi ~~ y/p) + \/V = Xi + y/p and 
X2 = X 2 , and dropping the unimportant constant term, we find 

X 2 ) = ixj 5 + ^>X\ - X 1 Xi . (178) 

The stability matrix $(Xi,X 2 ) at the critical point (X±,X 2 ) = (0,0) is 

(179) 



(0,0) 




Thus, at (xi,X2) = (^,0), as is mentioned above, the eigenvalue in the xi direction 
A X1 = 2y/p and the eigenvalue in the X2 direction A X2 = 0. For p > p , the two saddles 
and the local minimum are combined to form a single saddle at (xi, X2) = (y/P — o, 0), with 
(A X1 ,A X2 ) = 2(y/p — a, y/p — y/p — a). The X2 eigenvalue of the critical point (xi,X2) = 
(y/p — a, 0) changes sign as p passes through p . This implies a dual cusp catastrophe 
A_ 3 at (xi,X2) = (y/p, 0), where (p, z) = (po,0). To prove the existence of the cusp at 
(X11X2) = (y/p,®) rigorously, we need to reduce the catastrophe function Eq. 178 into the 
canonical form of the dual cusp A_ 3 : y/pX[ 2 + C 2 X 2 . Making a smooth change of variables, 
the same as Eq. 173, but, with x'% replaced by X[, and Xi replaced by JQ, we obtain the 
equation that the catastrophe function at (Xi,X 2 ) = (0,0) has to satisfy 

y/p-Xl - X X X\ + l -X\ = VpXf + C 2 X' 2 A = y/p-(Xl + 2X,A, + Aj) 

+C 2 (X 2 4 + 4X 2 3 A 2 + 6X^A 2 2 + 4X 2 A 3 2 + A 4 2 ) . (180) 

The equality is already satisfied for the terms of degree two. For the terms of degree three, 
we have 

-X X X 2 2 + X -X\ = 2 X /^X 1 A 1 = 2y/p-(A 20 Xf + A XX X\X 2 + A^X X X\) , (181) 
whose solution is 

A 20 = -— , A u = 0, A Q2 = ^-. (182) 
6y/p 2y/p 

For all the terms of degree four, we must have 

= 2 y /pX 1 k 1 + y/p-Al + C 2 X\ = 2y/p-(A 30 Xl + A 21 X*X 2 + A X2 X\X\ + A^X X X\) 

' -Xf - -^-=.X\X\ + -^-X\ + C 2 Xl , (183) 



36VP 1 &VP ^y/p 
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which implies 

A 30 = --!-, A 21 =0, A 12 = -±-, A 03 = 0, C 2 = --L. (184) 
72p Yip 4y/p 

Similarly, it can be shown that, Eq. 180 can be satisfied for the terms of all degrees greater 
than four. Thus, we have 

$ = Jp-X? - -±-X? . (185) 
A^p 

By rescaling the control parameters p^X[ = x' and 2~^p~^X' 2 = y', we obtain the canonical 
form of A_3 dual cusp catastrophe: $ = x' 2 — y' 4 . This proves the existence of the dual 
cusp at (Xi,X 2 ) = (0,0), or i.e. at (xi,X2) — (y/p,ty, an d hence, at p = p on the p-axis. 
For p > p , a second saddle also exists along xi < axis, at (xuX2) — (~VP ~ a ' 0) w ^ n 

( A xi> A x 2 ) = 2(- v / P- a >v / ^+v / P- a ); see Fi g- 7 - 

Although the study has been carried out on the symmetry axis z = to make the 

calculation simple, the procedure can be generalized to an arbitrary direction and the results 

are generic. When the point in control parameter space moves through the boundaries of the 

triangular region in an arbitrary direction, a fold or a dual cusp catastrophe occurs depending 

on whether the point moves through one of the three arc like sides or one of the three corners 

(dual cusps). Upon moving through a side line (fold) a minimum collides and annihilates 

one of the three saddles. At a fold, one of the three saddles and the local minimum are on 

top of each other. The saddle that is annihilated depends on which boundary side line the 

control point moves through. Upon moving through a dual cusp, two of the saddles and the 

local minimum collide, one of the saddles and the minimum annihilate each other, and one 

single saddle remains. At a dual cusp, two of the three saddles and the local minimum are 

on top of each other. Hence, the results are (qualitatively) independent of the restriction we 

made by focusing the p-axis only. The CDM flow described by the $ (Eq. 154) necessarily 

generates the tricusp caustic in physical space. In this sense the caustic rings are structurally 

stable. 



VII. CONCLUSIONS 

The infall of Cold Dark Matter (CDM) onto isolated galaxies such as our own Milky Way 
produces discrete number of flows and caustics in the halo CDM distribution. Caustics are 
locations in space where the density diverges in the limit of zero velocity dispersion. The 
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generic caustic is a surface at the boundary between two regions, one of which has n flows 
and the other n + 2 flows. In the limit of zero velocity dispersion, the density diverges as 
where a is the distance to the surface, on the side with n + 2 flows. If the velocity dispersion 
is small, but nonzero, this divergence is cut off, because, the location of the caustic gets 
smeared out. So the density at the caustic is no longer infinite, but merely very large. 
There are two types of caustics in the galactic halos: outer and inner. An outer caustic is 
a simple fold (A 2 ) catastrophe, located on a topological sphere enveloping the galaxy. An 
inner caustic is a ring whose cross-section is an elliptic umbilic (-D-4) catastrophe with three 
dual cusps (A_ 3 catastrophes); see Fig. 1. We call this cross-section a tricusp. Our focus 
was upon the caustic rings in this paper. In Sect. II, we reparameterized the CDM flow 
near a caustic ring, and obtained the ring equations in space, as single valued functions p(ip) 
and z(ip) of an angular variable ip parameterizing the tricusp. The density profile near the 
surface of a caustic ring is d(ip,a) = ^^-Q(a), where A(tp) is called the fold coefficient. In 
Sect. Ill, we derived the fold coefficient, everywhere on the caustic ring. A(ip) is minimum 
at the middle locations between the cusps. It increases monotonically as one approaches to 
the cusps, where it diverges in the limit of zero velocity dispersion. A caustic ring also has 
a specific geometry. Section IV is devoted to the differential geometry of the caustic ring 
surface. We obtained Gaussian, mean and principal curvatures as functions of ip. In Sect. 
V, we derived the gravitational lensing effects of the caustic rings for the line of sights that 
are parallel to the galactic plane of the ring and near tangent to the surface at any given tp. 
We estimated the image magnifications due to the caustic rings at cosmological distances 
and due to the nearby fifth ring of our own galaxy. For a given critical surface density, the 
effects are proportional to the fold coefficient A(ip) at the point where the line of sight is 
near tangent to the surface, and inversely proportional to the square root of the curvature of 
the surface along the direction associated with the line of sight. The magnification increases 
as the line of sight approaches to the cusps where it diverges in the limit of zero velocity 
dispersion. If the velocity dispersion is small but nonzero the cusps are smoothed out, hence 
the lensing effects are no longer infinite but merely large at the cusps. We used the lower 
and upper bounds of the effective velocity dispersions of the axion and WIMP flows in 
galactic halos to constrain the lensing effects at the cusps. For a cosmological axion caustic 
ring, we found that the magnification may range between 3% and 2800% at the outer cusp, 
and between 2% and 46% at the non-planer cusps. For a cosmological WIMP caustic ring, 
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on the other hand, we constrained the magnification between 3% and 28% at the outer 
cusp, and between 2% and 5% at the non-planer cusps. Because the upper bounds for the 
magnification at the cusps are obtained considering the minimum primordial value of the 
velocity dispersions of the CDM candidates in space, and because the observer's line of sight 
may not exactly be parallel to the galactic plane of the caustic ring in general, they should 
be regarded with precaution. In the limit of zero velocity dispersion, at a sample point near 
the outer cusps of cosmological CDM caustic rings, we found about 37% magnification. This 
effect is about 250 times greater than the effect obtained for the concave case considering the 
outer caustics, and about 37 times greater than the largest effect predicted considering the 
caustic rings where the case was convex, before. The nearby caustic ring, the fifth ring of the 
Milky Way, is only 1 kpc away (in the direction of observation) from us. We also estimated 
the magnifications of this ring of our own galaxy. Unfortunately, even near the outer cusp, 
the lensing effects are too weak to be observed with present instruments. In Sec. VI, we 
presented the correspondence of our formulation with the Catastrophe Theory. We derived 
the Catastrophe Function of the triaxial caustic rings and obtained the flow equations as 
the equilibrium points of this Catastrophe Function. The analysis of the Stability (Hessian) 
Matrix showed that the caustic rings are structurally stable; see also the recent simulations 
[20] confirming that the caustic ring is stable under perturbations. 



In this appendix, we show geometrically that the principal curvature radii on the caustic 
ring surface are the inverses of the corresponding principal curvatures. On the (x, y) and 
(y, Z)-planes of Fig. 8, the caustic cross-sections around an arbitrary point ip* on the surface, 
can be approximated locally by the circles, as indicated by the dashed arcs of the figure. 
On the (x, y)-plane, the circle that is tangent to the ring surface at the point (p(i[>*), z(ip*)), 
satisfies the equation 



where (p , Zq) is the center of the circle. By taking the first and second derivatives of both 
sides of Eq. 186 with respect to p, we obtain 



VIII. APPENDIX 



(p - p ) 2 + (z- z ) 2 = R\ 



'1p 1 



(186) 




(187) 
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FIG. 8: Same as Fig. 2 except that, now, we define a new Cartesian coordinate system (x, y, Z) 
at the arbitrary point ip*, such that (x,y) = (r), —ft). In the neighborhood of ip*, the caustic cross- 
sections in the (x, y) and (y, Z) planes can be approximated by the circles indicated by the dashed 
arcs. The circle in the (x,y)-plane has radius whereas, the one in the (x,Z)-plane has radius 



z- z = - 



'dV 
dp 2 , 



-i 



' dz" 
dp, 



(188) 



respectively. Using Eqs. 187-188 in Eq. 186, we find 



R 



1p 



± 



d 2 z 



dp 2 



-i 



1 + 



dz 
dp, 



(189) 



(p(V>*)> 2 W'*)) 

We may express the first and second derivatives of z with respect to p as 

dz z' d 2 z p'z" — z'p" 



dp p'' dp 2 



y3 



(190) 



where prime denotes derivative with respect to ip. Inserting, Eq. 190 into Eq. 189 and using 
the fact that p'z" — z'p" < on the whole surface, we obtain 

{ p ' 2 {^) + z' 2 {^))i l 



(191) 



z'(ip # )f/'(ip # ) - p , (ip*)z"(i/j„) /%(?/>*) ' 

In writing the last equality we compared with Eq. 112. 

We, next, want to show that R<p, the principal curvature radius in the plane defined by 
y and z = x x y (Fig. 8), is the inverse of k^. Around the origin of the (x,y, z)-coordinate 
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system of Fig. 8, which is chosen to be the point ip*, the caustic surface equation can be 
expressed as: 



y 



X 2 z 2 



+ 



2Rjf, 2R ( / ) 



(192) 



where > for < tp* < 2tt, whereas, R^ > only in the region < V** < 3 ■ -R^ < in 
the regions < ip* < and 4^ < ip* < 2ir; see Fig. 8. The intersection of the caustic surface 
with z = z*-plane is a circle of radius p*; therefore, near (p, z) = (p*,2*), the equation for 
this circle can be written as 

7.2 

(193) 



P~ P* = 



2p* 



On the other hand, from Eq. 40, in the z = z^-plane, we have 

x — (p - p*) cos#*, y — (p - p*) sin^„ , 



(194) 



where x — r], and y = — er. Inserting Eqs. 194 into Eq. 192, in the neighborhood of (p*, z*) 
we find that the intersection of the caustic surface with the z = ^-plane satisfies 

ry2 



y = (p - p*) sin 6** = 



2-fL, 



(195) 



where we neglected the term proportional to (p — p*) 2 . Combining Eqs. 193 and 195 we 
obtain, 

R<h 1 



sin 0* 

By the definition of the angle 9*, 

(p - p*) tanfl* = — (z - z*) , 



see Fig. 8. Hence, we can express tan#* = —^7, thus 



(196) 



(197) 



sin0* = ± 



[-091 


2 







(198) 



Insertion of Eq. 198 into Eq. 196 yields 

pCVQ^ty.) + **(&) j_ 

In writing the last equality we compared Ra, with Eq. 114. 



(199) 
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